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Abstract 

The Hubble Space Telescope (HST) , in its thirteen years of operation, has allowed us 
to observe properties of the universe humans have been able, until very recently, to 
probe only with their thoughts. This review presents a brief summary of a few of the 
highlights of HST discoveries, discusses their physical implications, and identifies 
unsolved problems. A broad range of topics is covered, from our own solar system 
to cosmology. The topics fall into the general categories of: planets (including both 
in the solar system and extrasolar), stellar evolution, black holes (including both of 
stellar-mass and supermassive) , galaxy formation and evolution, the determination 
of cosmological parameters, and the nature of the recently discovered "dark energy." 
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I. INTRODUCTION 



When Galileo Galilei peered through his small telescope 400 years ago, it resulted in an 
unprecedented period of astronomical discovery. There is no doubt that the Hubble Space 
Telescope (HST; Fig. 1) in its first 13 years of operation has had a similarly profound impact 
on astronomical research. But the Hubble Space Telescope did much more that that. It 
literally brought a glimpse of the wonders of the universe into millions of homes worldwide, 
thereby inspiring an unprecedented public curiosity and interest in science. 

The Hubble Space Telescope has seen farther and sharper than any optical/UV/IR tele- 
scope before it. Unlike astronomical experiments that were dedicated to a single, very 
specific goal (like the Cosmic Background Explorer or the Wilkinson Microwave Anisotropy 
Probe), the Hubble Space Telescope's achievements are generally not of the type of sin- 
gular discoveries. More often, the Hubble Space Telescope has taken what were existing 
hints and suspicions from ground-based observations and has turned them into certainty. In 
other cases, the level of detail that HST has provided forced theorists to re-think previous 
broad-brush models and to construct new ones that would be consistent with the superior 
emerging data. In a few instances, the availability of HST's razor-sharp resoution at critical 
events provided unique insights into inividual phenomena. In total, by observing tens of 
thousands of astronomical targets, the Hubble Space Telescope has contributed significantly 
to essentially all the topics of current astronomical resarch, covering objects from our own 
solar system to the most distant galaxies. 

The 12.5-ton orbiting observatory was launched into orbit on April 24, 1990. It circles the 
Earth every 90 minutes in a 330 nmi (607 km) orbit and operates around the clock above all 
but the thinnest remnants of the Earth's atmosphere. The telescope provides information 
to many individual astronomers and teams of scientists worldwide, and it is engaged in the 
study of virtually all the astronomical constituents of the universe. 

Crucial to fulfilling the objective of a 20-year mission is a series of on-orbit manned ser- 
vicing missions. The first servicing mission (SMI) took place in 1993, the second (SM2) 
was flown in 1997, and the third was separated into two parts, one (SM3A) in 1999 and the 
second (SM3B) in 2002. A fourth mission is currently planned for the end of 2004 (pend- 
ing the investigation into the shuttle Columbia disaster). During these missions, shuttle 
astronauts upgrade the the observatory's capabilities (by installing new instruments) and 
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perform planned maintenance activities and necessary repairs. To facilitate this process, 
the telescope was designed so that its science instruments and several vital engineering sub- 
systems have been configured as modular packages with standardized fittings accessible to 
astronauts in pressurized suits. 

Hubble was designed to provide three basic capabilities, all of which are essential for 
innovative astronomical research: 

1. High angular resolution — the ability to image fine details. The spatial resolution of 
HST in the optical regime is about 0.05 arcseconds. 

2. Ultraviolet performance — using the fact that the observatory is above the Earth's 
atmosphere to produce ultraviolet images and spectra. The Hubble Space Telescope is 
sensitive down to a wavelength of 1150 A (below the Lyman Alpha line of hydrogen). 

3. High sensitivity — the ability to detect and take spectra of even very faint objects. 
The Hubble Space Telescope is about a hundred to a thousand times more sensitive 
in the ultraviolet than the previous space observatory (the International Ultraviolet 
Explorer) . 

Astronomers and astrophysicists using HST data have published over 3000 scientific papers 
to date. Hubble's findings are thus far too numerous to be described even briefly in one 
article. The following sections therefore simply highlight a few of the astronomical discov- 
eries which, in my clearly biased view, have significantly advanced our understanding of the 
cosmos. In the spirit of progressing by "powers of ten" in cosmic distance, I shall start with 
the solar system and finish with cosmological distances. I apologize in advance to the many 
astronomers whose important work will not be presented in this article. I also regretfully 
acknowledge the incompleteness of the list of references. My list should be regarded as rep- 
resentative rather than comprehensive. Finally, the different sections are not intended to be 
exhaustive reviews — they are concise summaries of HST contributions to major astrophysical 
poblems. 

II. OUR OWN BACKYARD THE SOLAR SYSTEM 

The Hubble Space Telescope has proved well suited for the study of objects in our own 
solar system. In fact, it has produced images of the outer planets that approach the clarity of 
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those obtained from Voyager fiybys. The Hubble Space Telescope's images of Mars have also 
been surpassed only by close-up shots taken by visiting space probes. The Hubble Space 
Telescope's images of the auroras of Jupiter and Saturn, of the disintegration of comet 
LINEAR, of the binary Kuiper belt object KBO 1998 WW31, and of a global dust storm on 
Mars have been truly spectacular and informative. However, no observation could compete 
with the drama provided by Comet Shoemaker-Levy 9 (SL-9) slamming into Jupiter. 



A. The collision of Comet SL-9 with Jupiter 



From a cosmic perspective, the impact of the ninth periodic comet discovered by the 
team of Carolyn and Gene Shoemaker and David Levy on Jupiter was unremarkable. The 
cratered surfaces of many planetary bodies and satellites and the prehistoric mass extinctions 
on Earth are a testament to the fact that such impacts are ubiquitous in the history of the 
solar system. Furthermore, even the "seeds" of life on Earth may have been planted by such 
a bombardment. 

From a human perspective, however, this collision marked a "once in a lifetime" event, 
which caused the mobilization of all the astronomical resources for an unprecedented obser- 
vational campaign in July 1994. 

The Hubble Space Telescope played a crucial role in the SL-9 campaign, both in the 
initial characterization of the comet fragments and in producing unique images and spectra 
of the impact events themselves and later of the impact sites. 



The multiple comet SL-9 was discovered in late- March 1993 (Shoemaker et all [1993) 



Subsequent images taken with HST in July of 1993 ((Weaver et al. 



19941 1 revealed that the 



comet had spectacularly split into 20-odd fragments, probably because of tidal disruption 
(during a previous close approach to Jupiter), producing a train of fragments fittingly dubbed 
a "string of pearls" (Fig. 2). 

The conditions for tidal break-up of a self-gravitating incompressible cometary sphere 



can generally be written in the form (e. g. . iDobrovolskis . 
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1990; 



Sekanina et al. 



19941 ) 



Pc ( d 
P P \Rp, 



< 
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ll 



(1) 



where p c , p p are the densities of the comet and planet, respectively, R p is the planet's radius, 
d is the distance of the comet from the planet, P c is the comet's central pressure, U is a 
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measure of the comet's tensile strength, and k is a constant of order unity. The facts that: 

, and (ii) the grad ual disappearance of 



(i) secondary fragmentation was observed by HS 
condensations was also well documented by HST ((Weaver et a/I Il995f ). argued strongly for 
tidal splitting of a discrete nucleus as a result of gradual fissure propogation. Computations 



of the orbit suggested that th e comet would collide with Jupiter in July 1994 (|Nakano 



Yeomans and Chodas 



1993). 



1993; 



The impact of the first fragment (fragment "A") occurred on July 16, 1994, followed by all 
the other fragments smashing into Jupiter's atmosphere over the following w e ek. T he Hub- 



ble Space Telescope's images taken during the impact events (jHammel et al 



199.4 



showed 



spectacular plumes resembling nuclear "mushrooms" above the limb of Jupiter (Fig. 3). 
Typically, the plumes attained altitudes of about 3000 km within 6-8 minutes of impact, 
and were seen falling and spreading within 10 minutes of impact. The manner in which 
the height of the plume scales with SL-9 parameters and the ambient pressure can be ob 



tained from the following simple c onsiderations (e.g. 



1995 



Zahnle and MacLow 



Chvba et al 



1993 



Field and Ferrara . 



1994). Once aerodynamic forces overcome material strength, 



the incoming fragment deforms to form a flattened "pancake" due to the difference in ram 
pressure acting on the front and on the sides. Due to the increase in the cross-section (and 
concomitantly in the drag), the fragment is brought to an abrupt halt, accompanied by an 
explosive energy release. If we denote the energy being released per steradian by Ei/Air, 
then the obtained ejection velocity is given by 



V ■ ~ 



E ig 



(2) 



/ A7rp a H 3 V 4vrP a # 2 

where p a , P a are the ambient density and pressure, respectively, H is the pressure scale height 
in the atmosphere and g is the gravitational acceleration. The plume height, therefore, scales 
like 

V e 2 Ei 

Vume - ~ 77T^- ( 3 ) 



Detai 



le.g. 



2g P a H* 

ed numerical s imulations have confirmed the basic physical picture described above 



MacLow 



1996). 



A combination of the observations from the Keck Telescope, HST, Palomar, and the 
Galileo spacecraft, eventually produced a co mprehensive explana tion for the observed 2.3yU 
lightcurve (shown in Fig. 4 for fragment R; iGraham et all Il995). The sequ ence of events 



can be described phenomenologically by the cartoon in Fig. 5 flZahnle , 



1996J) 



The first precursor marked the meteor trail (with a brightening timescale of order H/v ~ 
Is). As explained above, most of the fragment's kinetic energy was released at the last second 
of its existence, resulting in an explosive fireball. The second precurso r marks the ris e of th e 



fireball above Jupiter's limb, having already cooled to ~ 500-700 K ((Carlson et al. 



1995). 



As the fireball expanded, silicates condensed, followed by carbonaceous matter. The sunlit 
parts of these condensates were the tracers that marked the plume seen by HST. The main 
peak in the light cury e was caused by the ejecta plume falling back onto the atmosphere 



( Nichol son et al. 



1995). The rise and fallback occurred on a timescale of JSh p \ nme / g ~ 10 3 s 



(for a plume height of ~ 3000 km). The final, smaller peak, was probably caused by material 
bouncing off the atmosphere and falling back. 

Another interesting phenomenon revealed by the HST images of the impact sites was the 
propagation of waves (Fig. 6). The radius of the pro minent dark ring in F ig. 6 is 3700 km, 



and a smaller ring, of radius 1750 km is also visible (|Hammel et al 



1995L 



Measurements of the radii of the circular rings as a function of time showed that the "fast" 
wave was seen to spread at a constant speed of 450 m s _1 , while the speed of the "slow" wave 
(corre sponding to the fainter ring) was somewhat less constrained. Ilngersoll and Kanamori 
f)l996l ) examined many possibilities for the identification of the fast wave, including /-modes, 
p-modes, r-modes and g-modes. These waves are classified by the different restoring forces 
that cause them to propagate. In acoustic waves (p-modes) the restoring force is compress- 
ibility; in surface waves (like those in the ocean; /-modes) it is gravity; in inertial oscillations 
(r-modes) it is the Coriolis force, and in internal gravity waves (g-modes) it is gravity. Inger- 
soll and Kanamori showed that the only waves that match the observed speed of 450 m s _1 
are g-modes propagating in a water cloud which has an O/H ratio that is ten times higher 
than that ratio in the Sun, eH 2 o = 10 (the waves were too slow to be acoustic waves and 
moved at too constant a speed to be rings of debris). The propagation speed of these modes 
varies as the square root of the water abundance. The requirement of a high O /H ratio was 
found to be in conflict with data from the Galileo probe, which suggested that from above 
the cloud tops down to the 20 bar level the water was less abundant than in the Sun. The 
cause for this discrepancy is still unknown, although it has been suggested that the probe 
may have entered a clearing in the clouds (a "5-micron hot spot") associated with a dry 
downdraft. This uncertainty demonstrates that even the highest quality data do not provide 
all the answers. Nevertheless, there is no question that the impact of comet SL-9 on Jupiter 
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provided a unique opportunity to study such events. 



III. EXTRASOLAR PLANETS 



The solar system was not the only place to provide for planetary drama. During the past 
eight years we have witnessed the number of known extrasolar planets orbiting Sun-like stars 
going from zero to about one hundred at the time of this writing! The human fascination 
with the possibility of discovering extraterrestrial life has helped turn planet detection into 
one of the major frontiers of today's astronomy. 

The Hubble Space Telescope's contributions to this exciting field have been related mostly 
to the discovery of circumstellar disks from which planets form, and these will be described in 
Section IV. In two cases, however, HST has produced three truly unique sets of observations, 
one in a globular star cluster, and two following the light curve and spectrum of a transiting 
planet. 



A. Where have all the planets gone? 



The discovery of a planet orbiting the star 51 Peg (M ayor and Queloa . 1199 5) in a re- 
markably tight orbit of 4.231 days challenged prevailing theoretical views and marked the 
start of a "Golden Age" in planet discoveries by radial velocity surveys (see, for example, 



Marcv et al. 



2000). The existence of gas giant planets in very close orbits (at 0.04-0.05 AU) 



enables even photometric searches for these planets ("hot Jupiters"). At such small sepa- 
rations, the probability of the planet transiting ( "eclipsing" ) the parent star, given random 
orbital inclinations, is about 10%. More precisely, from Kepler's laws, the duration of a 
transit in hours (reduced by ir/4 for average chord length) is given by 



Ttr 



lAM'^R^ hr 



(4) 



where M*, R* are the stellar mass and radius in solar units and P or b is the orbital period in 
days. The probability of a transit per system is given by 



Pr 



tran 



~ 23.8%M- 1 / 3 J R*P"b /3 



(5) 



Transits are very important because of the following reasons: (a) With transits one can 
immediately obtain the size of the planet [the transit depth is equal to (i? p /-R*) 2 ]. (b) In 
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the case of transits, measurements of radial velocities immediately provide the planet's mass 
(and not just m p sin i, where i is the otherwise unknown orbital inclination), (c) A knowledge 
of the mass and radius can be used to test theoretical mass-radius relations for close-orbiting 
planets. 

In 1998, a team of astronomers led by Ron Gilliland of the Space Telescope Science 
Institute embarked on the ambitious program of searching for transits in some 40,000 stars in 
the core of the globular cluster 47 Tucanae (Fig. 7). This cluster is a particularly interesting 
target for such a search for several reasons, (i) For a typical star in the cluster, M* = 0.81, 
i?* = 0.92, and a typical "hot Jupiter" orbital period (e.g., P or b = 3.8 days), the probability 
of transit is 9.6% and the duration is 2.2 hours. For a planet of radius R p = 1.3-Rj (where 
R.j is Jupiter's radius, as observed for the transiting system HD 209458; 



2000; 



Henry et al 



Charbonneau et al. 



20001 ). the expected transit depth is about 2%, easily detectable by HST. 



(ii) Given that the frequency of "hot Jupiters" with orbital periods shorter than 5 days (ten 
are known at the time of writing) in the solar neighborhood is about 0.8%-l%, with a 10% 
chance of transit one expects roughly one in every 1000 surveyed stars in 47 Tuc to show a 
transit (and therefore a few tens of detections when monitoring ~ 40,000 stars), (iii) The 
globular cluster 47 Tuc is about 11 billion years old and it has a total met allicity (abundance 
of ele ments heavier than helium) of about one-third that of the Sun ((Salaris and Weiss . 



19981 ). The frequency of planets in such an environment can therefore have important 



implications for models of planet formation. 

The Hubble Space Telescope observed 34,091 stars in 47 Tuc, obtaining time series pho- 
tometry over a period of 8.3 days (continuously, except for Earth occultations and passages 
through the South Atlantic Anomaly). A more detailed calculation of the expected number 
of transi ts, taking into accoun t the actual distribution of V magnitudes of the stars in the 
sample (Gi lliland et all |2000J) and time series noise, showed that 17 transits should have 
been seen, assuming the same frequency of hot Jupiters as in the solar neighborhood. None 
was detected! Undoubtedly, this was not a problem with the detection method. As many 
as 75 variable stars have actually been detected. Furthermore, an eclipsing system with a 
period of 1.34 days and an eclipse depth of ~ 3% was actually discovered. However, the 
presence of a secondary eclipse in the system showed that this was a grazing eclipse by a 
low- mass (K dwarf) star and not a planet. 

The direct implication is therefore that with a very high degree of confidence the frequency 
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of hot Jupiters in 47 Tuc is at most 1/10 that in the solar neighborhood. 

The main question is then: why are there much fewer orbiting planets in 47 Tuc? There 
are, in principle, two possible explanations for the paucity of planets in orbits around globular 
cluster stars: (i) Planets do not form at all (or form much more rarely) in globular clusters, 
(ii) Planets are torn from their parent stars. 

Planet formation may be suppressed in old stellar clusters because of two main reasons: 
low metallicity, or photoevaporation of protoplanetary disks. 

Conventional wisdom of giant planet formation suggests that this is a two-step process. 



In the first, the collisional accumulation o 
away growth of a solid core (e.g. 



Lissauer 



' rock y and big planetesimals leads to a run- 
1987) of about ten Earth masses (10 M e ). In 



the second, the core acquires a gaseous atmosphere, which eventually collapses, leading to 
the rapid accretion of hydrogen and helium from the protoplanetary disk. One might sus- 
pect, therefore, that in low-metallicity environments, the nucleation of dust grains would be 
greatly suppressed, and concomitantly planets would form more rarely. Furthermore, it has 
generally been noted that s tars which bear planets are, on average, enriched in metals i n 



Laughlin 



2000; 



2002). 



1 1 1 

comparison to the Sun (e.g.. lButler et all 120001 : lGonzalea . ll997 
and this trend is even more pronounced for the short-period planets. Thus, the paucity of 
hot Jupiters in 47 Tuc may simply reflect the cluster's low metallicity. 

There is another reason why giant planet formation may be suppressed in dense stellar 
clusters. The disk around any star in a cluster is exposed to ultraviolet radiation from 
massive cluster stars. This radiation heats the disk surface, raises the sound speed, c s , and 
causes gas beyond a radius of 



Rg~ 



0.5- 



.GM, 



(6) 



to become unbound and flow a way as a thermal wind. The mass-loss rate in _ such a photo- 
evaporation flow scales as (e.g., 



Rertoldi and McKeel. |l990; 



Johnstone et al 



1998) 



outflow 



&l 2 d- 1 



(7) 



where $ is the flux of ionizing photons and d is the disk's distance from the emitting source. 
In a star forming region like Orion, the ionizing flux is of order $ ~ 10 49 s _1 , of the same 
order as the flux expected in a cluster with n c i us ter ~ 10 3 members. The mass loss rate from 
one of the disks in Orion, HST 182—413, has been estimated from HST observations to be 



^outflow ~ 4 x 10 7 M yr 1 . Since the mass of that disk is estimated to be ((Johnstone et al 
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19981 ) Mdisk ~ 0.04 M Q , the disk lifetime is expected to be Tdi S k ~ Afdi S k/M out fl ow ~ 10 5 yr. 
For somewhat smaller disk sizes (e.g., like that in the solar system ~ 30 AU) the disk lifetime 
can be somewhat extended (the EUV-induced flow scales as r d ( sk ). Nevertheless, the disk 



lifetime may be significa ntly shorter t 
order 10 6 -10 7 yrs (e.g. 



r an th e planet formatio n timescale, e stima ted to be of 



Pollack et al 



1996 



Shu et al 



1993J). lArmitagel (fcOQOh calculated 



the expected disk lifetime as a function of the cluster richness. He found that giant planet 
formation is strongly suppressed in custers with ~ 10 5 stars out to around 1 pc or more. 

There are two ways in which, in principle at least, the conclusion of a low probability 
for planet formation in clusters, due to disk destruction, may be avoided: (i) If giant planet 
formation occurs not via a core accretion mechanism, but rather through direct and rapid 
(within ~ 1 3 yr) fragmentation (due to a gr avational instability) into self-gravitating clumps 



of gas (e.g., iBoss 



2000; 



Maver et al 



20021 ) . (ii) If the formation of high- mass stars (which 



produce the photoevaporating radiation) is significantly delayed (by ~ 10 6 yr) compared to 
the formation of low mass stars. At present both of these processes are sufficiently uncertain 
that no definitive conclusion about their viability can be drawn. 

The second possible explanation for the absence of "hot Jupiters" in 47 Tuc is that even 
if giant planets do form, they are torn away from their parent stars because of encounters 
in the crowded cluster environmen t, and thus cannot exh i bit tra nsits. This possibility has 



been investigated in some detail by 



Davies and Sigurdssonl (j200fl ) . The cross-section for two 



stars with a relative velocity at infinity to pass within -R m i n from each other is given by 

<7 = 7ri^(l + ^J , (8) 



where V is the relative velocity (V 2 oc (Mi + M 2 )/R Ta i n ) at closest approach. The second 
term expresses the effect of gravitational focusing. When V ^> V^, as is expected in globular 
clusters, a oc -R m i n . The timescale for an encounter between a star-planet system and a star 
is of order (r cnc = 1/naV) 



r enc ~ 3 x 10 8 yr 



a 



d 



M 



(9) 



v io 5 P c- 3 ; vio km s"V V 1 au; V M 

where n is the number density of stars of mass M and d is the semimajor axis of the planetary 
orbit. Encounters therefore occur on a timescale that is much shorter than a Hubble time 
in a cluster like 47 Tuc, which has a density of about 1.5 x 10 5 pc~ 3 at its cent er and about 
an order of magnitude lower at its half-mass radius (e.g., iHowell et all I2000T ) . Generally, 
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binaries (with masses M\ and M 2 for the primary and secondary, respectively) in which 
the binding energy is lower than the kinetic energy of the colliding star (of mass M3) are 
disrupted. The condition for breakup is 

2 > GM 1 M 2 (M l + M 2 + M 3 ) 

M 3 (M 1 + M 2 )d ' 1 J 

Since M 2 <C Mi,M 3 , in the case of a star-planet system, breakup of the system is ex- 
pected to occur with a high probability when the third star passes within a distance d. 



Davies and Sigurdssonl (|200lh performed a large number of simulations for various values 



of Ko/Vorb and considered both encounters (of star-planet systems) with single stars and 
with binaries. They found that wide (rf^0.3 AU) planetary systems are likely to be broken 
up within the half-mass radius of 47 Tuc, but tighter systems (d^O.l AU) or systems in 
less dense regions may survive. The results of Davies and Sigurdsson still do not provide a 
complete explanation for the absence of transits in 47 Tuc. Future observations, perhaps of 
less dense, somewhat higher metallicity clusters (like NGC 6352) will be required to explain 
which factor is responsible for the dearth of planetary transits in 47 Tuc. An attempt to 
detect transits associated with Galactic bulge and disk stars (that span 1.5 dex in metallic- 
ity, and are in an environment that is orders of magnitude less dense than the center of a 
globular cluster) would also be very valuable. An understanding of the environments that 
are either conducive to or prohibitive for the existence of planetary transits may constitute 
an important step in the study of planet formation. 



B. The transiting planet HD 209458 



The low-mass companion to the star HD 209458 was the first extrasolar planet 



to transit the optic al disk of its parent star ()Charbonneau et al 



Mazeh et al. 



2000; 



Henry et al 



bund 



2000; 



2000). As I explained in the previous section, transits offer a unique opportu- 



nity to deter mine the properties o f the orbiting planet. The Hubble Space Telescope followed 



four transits iBrown et al. 



2001) and obtained an exquisitely detailed light curve (Fig. 8; 



note that the depth of the eclipse is only 1.7%!). The orbital period was determined to be 
3.5247 days, the lower limit on the mass of the planet M p sini = 0.69 ± 0.05M/ (with the 
inclination angle % = 86°6±9°14), and the radius of the planet R p = 1.347 ±0.060 Rj. Most 
impressively, the precision of the HST light curve allowed even for searches for rings around 
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the planet, for constraints on planetary satellites, and for probing the planet's atmosphere. 

A ring system with significant opacity around the planet would cause distortions of the 
light curve relative to that of a spherical body. In particular, one would expect small dips in 
the light curve before the first and after the fourth contact. The observations were consistent 
with no rings (maximum ring radius consistent with observations was 1.8 R p ). 

A satellite orbiting the planet might be detectable either from its photometric signature 
(the satellite would block light in addition to t hat obstructed by t he planet), or from its 
influence on the orbital motion of the planet. iBrown et al\ (|200l[ ) showed that satellites 
larger than 1.2 (where R© is the Earth radius), or with masses larger than 3 M e are 
excluded by the data. 

Finally, precision spectrophotometric observations in the region of the sodium reso- 
nance doublet at 589.3 nm revealed that the photometric dimming during transit was 
deeper by (2.32 ± . 57) x 10~ 4 relative to simultaneous observations in adjacent bands 



(Ch arbonneau et al. 



2002). This additional dimming has been interpreted as absorption by 



sodium in the planet's atmosphere. In fact, the existence of a detectable sodium featur e 
in the spectrum has been predicted by theoretical models (e.g.. Iseager and Sasselov . 2000 ). 
Furthermore, observations of brown dwarfs (stellar objects below the hydrogen burning 
limit) with similar su rface temperatures show strong absorption in alkali metal lines (e.g., 



Burrows et al 



20001 ) . HD 209458 therefore represents the first detection of an atmosphere 



of an extrasolar planet. While not a shocking result in itself, since giant planets are clearly 
expected to have atmospheres, the actual detection marks the beginning of a new era in 
extrasolar planets research. 

Even more interestingly perhaps, HST observations of HD 209458 during three tran- 
sits revealed atomi c hyd rogen absorption of ~15% over the stellar Lyman a line 
( Vidal-Madiar et all |2003). A comparison of the observations with models showed that 
the observations could be explained in terms of hydrogen atoms evaporating and escaping 
the planet in an asymmetric cometary-like tail. To account for the observed absorption 
depth, the simulations implied a minimum escape flow rate of ~ 10 10 g s _1 . 

In the future, other spectral features (e.g., water vapor; with implications for the searches 
for extraterrestrial life) can be searched for. Furthermore, in principle, HST can detect even 
the secondary eclipse (when the planet is eclipsed by the star), expected to be characterized 
by a ~ 10 -5 decrease in the luminosity. This will allow for a determination of the planet's 
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albedo, and for a comparison with models of planetary atmospheres. All of these results, 
whether already obtained, or expected, are somewhat of a surprise, since HST has been 
heralded primarily as a "cosmology machine," not as an exquisite tool in exo-planetary 
research. 

Moving now from planets to stars, the Hubble Space Telescope has been particularly 
instrumental in revealing new details of stellar births and deaths. 



IV. STARS AT BIRTH AND DEATH 



Using HST's high resolution on one hand and its infrared capabilities on the other, as- 
tronomers have been able to probe the dusty environments of young stellar objects (YSOs) 
and of stars in their late evolutionary stages with unprecedented detail. While these ob- 
servations have not revolutionized the field, they have, on one hand, led to a much deeper 
understanding of the processes involved, and on the other, opened an entirely new set of 
questions, while producing some of the most spectacular images. 



A. Outflows and jets from young stellar objects 



Molecular clouds are the reservoirs of mass and angular momentum from which stars 
are born. In the initial phase, dynamical infall occurs. The finite angular momentum of 
infalling cloud material leads to the formation of an accretion disk. Angular momentum 
is transported outward in this disk resul ting in the accret i on of mass (and some angula r 
momentum) onto the central object (e.g., 



Camenzinc 



mdl . 199(1 



Naiita . 



2000; 



Shu et al. 



1994). 



For a long time it has been known that outflows and collimated jets are signposts of stellar 



Konigl, 



198 



Reipurth and Heathcote 



19931 ). In particular, many 



birth (see, for example, 
of the radiative shocks known as Herbig-Haro (HH) objects were f ound to be assoc iated with 
highly collimated jets emanating from the vicinity of YSOs (e.g., 



Reiourth 



cases, these jets were found to be several parsecs in length (jReipurth et al 
shows a few remarkable HST images of jets from YSOs. 



1999). In some 
Figure 9 



In spite of the ubiquity of astrophysical jets, we s till lack a co mprehensive theory of 



their acceleration and collimation (see e.g., review bv iLivio . 2000). The most promising 
mechanism relies on an accretion disk (around a central compact object like a YSO or a black 
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hole) threaded by a large-scale poloidal magnetic field. Some magnetic flux is assumed to be 
in open field lines, inclined by an angle of more than 30° to the vertical to the disk's surface. 
Ionized gas is forced to flow along field lines. Since the foot points of these lines are anchored 
in the disk and rotate with it, material is accelerated cen trifugally like beads on rotating 



wires (jBlandford and Pavne . 



1982 



Qgilvie and Livid . 1200 II ). In this picture, the acceleration 



basically stops at the Alfven surface, where the kinetic energy density in the jet becomes 
comparable to the magnetic energy density. Collimation, however, occurs primarily outside 
the Alfven surface. There, the field gets wound up by the rotation, generating toroidal 
loops. The curvature force exerted by the toroidal field acts in the direction of the rotation 
axis, thus achieving collimation by "hoop stresses." Alternatively, in a vertical field of the 
form B z ~ (r 2 /R- m + l) -1 / 2 (where R m is the radius at the accretion disk's inner edge), in 
whic h the flux is la rgest at the outer disk, the field lines have a naturally collimating shape 



(.e.g. 



Spruid . I1996T ). Physically, poloidal collimation is achieved in this the materia. 



encounters the high flux in the outer disk (especially if -Rdisk ~ -^Alfven)- 

While this theoretical picture of jets being accelerated and collimated by accretion disks 
has been largely developed prior to any HST observations, HST has provided the first direct 
evidence for the fact that jets indeed origin ate at the centers of d isks (at least in the case of 
YSOs). Observations of the objects HH 30 ((Burrows et allied . Fig. 10), DG Tau (Fig. 9), 
and a few other YSOs clearly show the jet emanating from the disk center, and they reveal 
even the illuminated upper and lower disk surfaces. 

Furthermore, the high resolution of HST has enabled the determination of proper m otions 
in YSO jets, from images obtained over a time interval of a few years (e. g. . iBallvl 120031 ) . The 
highest velocities have been observed along the jet axis. Typically, the velocities are of order 
200 to 400 km s -1 . Since these are precisely of the order of the Keplerian velocity in the inner 
disk around a YSO, the proper motion observations provide additional support for the jet 
formation scenario described above. Incidentally, similar observations of the optical jet in the 
active galaxy M87 also showed proper motion. In this case at the apparent "superluminal" 
speed was of 4c-6c ( Biretta. Sparks, and Macchettol . 1999 ). again corresponding to the fact 
that one expects relativistic speeds from the vicinity of a central black hole. 

Another related phenomenon th at only HST could discover is the existe nce of externally 
irradiated small jets ("micro jets' - ' 



Ballv. O'Dell and McCaughrean . 



2000). Many of these 



microjets (Fig. 11) are only about 0.1" wide, and are therefore usually undetectable against 
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the nebular background in ground-based observations. For irradiated jets, observations of 
the Ha surface brightness, the emission measure (EM = J n 2 e dx\ where n e is the electron 
density and x is the linear size of the emitting region), and the jet width, allow for a 
determination of the electron density. Together with the jet speed one can therefore obtain 
the rate of mass loss in the jet, Mj. Typically, the irradiated small jets in the Orion nebula 
are characterized by Mj ~ 10~ 9 M Q yr _1 , at least an order of magnitude lower than the 
rates observed in the long jets associated with Herbig-Haro objects. Generally, the rate of 
mass loss in jets from YSOs is found to be about 1-10% of the rate at which mass is accreted 
through the disk onto the young stellar object. 

In addition to the important disk-jet connection, HST observations of YSOs have provided 
another interesting new element, this time related to planet formation. 



B. Protoplanetary disks 



Contrary perhaps to the expectation that protoplanetary disks would be deeply embed- 
ded within the clouds from which they form, and they would therefore be inaccessible to 
optical observations, HST revealed many dozens of protoplanetary disks ("proplyds"; e.g 



Ball v. O'Dell and McCaughrean 



2000; 



O'Dell 



2001 



O'Dell et al, 



1993 



O'Dell and Wong . 



19961 following the initial correct identification by Churchwell et al. 1987 and Meaburn 
1988). Many of these disks are seen silhouetted against the background nebular light (when 
they are shielded from photoion i zation) , with some possessing io nized skins and tails (e.g., 



Ballv. O'Dell and McCaughrean 



2000; 



Hennev and O'D ell. 



19991 . Fig. 12). 



The ubiquity of the protoplanetary dust disks (they are seen in 55%-97% of stars; 



Hillebrand et al 



1998; 



Lada et al. 



2000) demonstrates that at least the raw materials for 



planet formation are in place around many young stars. Indeed, in a few cases, like the 
dust ring and disk in HR 4796A and the nearly edge-on disk surrounding Beta Pictoris, the 
detailed HST image s reveal gaps and warping (respectively) t hat could represent the effects 



of orbiting planets fK alas et al 



2000; 



Schneider et al. 



1999). 



Another aspect of the protoplanetary disks that is significant for planet formation is the 
discovery of evaporating disks in the Orion Nebula. As was noted in Section IIIA, some of 
the Orion proplyds were shown to be evaporating (due to photo-ablati on by UV radiation 
from young, nearby stars) at rates of M ~ 10~ 7 to 10 -6 M yr _1 (e.g., iHennev and O'Dell . 
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19991 ) . Given that the masses of these disks are typically of order 10 2 M (if normal 



interstellar grains are assumed, so that the observed dust emission can be scaled to the total 
mass), this implies lifetimes for these disks of 10 5 years or less. There exists, however, some 
evidence that the grain s izes in Orion's disks may, in fact, be relatively large — perhaps of 
the order of millimeters ((Throopl . |2000J). The latter conclusion is based on the fact that 
the outer portions of the disks appear to be gray (they do not redden background light), 
and on the failure to detect the disks at radio wavelengths in spite of the implied large 
extinction in the infrared (hiding the central star in some cases) . The observations are thus 
consistent with grain sizes in excess of the radio wavelength used, of 1.3 mm. When we think 
about the potential implications of these two findings (about disk lifetimes and grain sizes), 
we realize that they may have interesting consequences for the demographics of planets 
in Orion. The relatively short disk lifetimes but relative large grain sizes may mean that 
while rocky (terrestrial) planets can form in these strongly irradiated environments, giant 
planets (that require the accretion of hydrogen and heliu m from the protoplanetary disk ) 
cannot (unless their formation process is extremely fast; iBossl 120001 : iMaver et all |2J302) . 
It is nevertheless clear from the many observations of "hot Jupiters" (giant planets with 
orbital radii ;$ 0.05 AU) that less extreme environments do exist, in which giant planets not 
only form, but also have sufficient time to gravitationally interact with their parent disk and 
migrate i nward, to produce the distribution in orbital separations we obse rve today (see, for 



example, 



Armitage et a,l 



2002; 



Lin. Bodenheimer. and Richardson . 



1996). 



While disks around young stars produce jets and form planets, similar structures around 
old stars help perhaps to shape incredible "sculptures" around dying stars. 



C. Morphology of stellar deaths 

The evolution of stars is determined primarily by their mass, because evolutionary 
timescales are set by the rate of consumption of the nuclear fuel. Figure 13 shows the 
evolutionary tracks of 1 M Q , 5 M , and 25 M stellar models in the luminosity-effective 
temperature plane. It is clear that certain phases of stellar evolution are characterized by 
considerable mass ejection, either via relatively gentle processes such as radiation pressure 
on atomic lines or on dust, or through stellar explosions. A process of the former type pro- 
duces the objects known as planetary nebulae — shells of ejected matter from giant stars that 
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are ionized by the radiation emitted by the hot (^,30,000 K) central core. The explosions 
known as Type II supernovae are the result of the dynamical collapse of the iron cores of 
massive stars, and the ejected nebulae form the objects known as supernova remnants. 

One of the most striking phenomena revealed by HST observations of stellar deaths and 
of stars in their late stages of evolution is the fact that axisymmetric nebulae are extremely 
common. This observation remains true even after recognizing that the morphologies are 
blurred by projection effects, and by the lines used to obtain the image (i.e., the [O III] 
image may look different than the [N II] image). Furthermore, axisymmetry is found not 
only in planetary nebulae (PNe), that are formed by intermediate-mass stars (1-8 M ), but 
also around Luminous Blue Variables (LBVs), that are extremely massive stars (~ 100 M Q ), 
and supernovae (SNe), that represent the death of massive stars (8-30 M ). 

Eta Carinae (rj Car), for example, belongs to the small class of very massive stars known 
as Luminous Blue Variables. These stars are believed to represent a rapid evolutionary 
phase, in which the stars experience severe mass loss (losing many s o lar m asses in ~ 10 4 



years), sometimes via giant outbursts (e.g., 



Davidson and Humphriesl . |1 



In the 1840s 



7] Car suffered such an outburst, increasing in brightness by several magnitudes and ejecting 
a large amount of material. The star stabilized around 1870 (except for a minor eruption 
between 1887 and 1895). The precise cause for this giant outburst is still unknown, but it has 
produced a spectacular, bipolar nebula commonly referred to as the Homunculus (Fig. 14). 
In addition to the Homunculus, the HST images revealed for the first time the presence 
of a ragged ejecta disk around the "waist" of the hourglass structure, composed, at least 
partly, of radial streaks. The total mass in the equatorial disk has been estimated (based 
on conventional gas to dust ratios) to be at least 0.1-0.2 M , but this could underestimate 
the mass significantly. 

Recent spectroscopic observations of rj Car by the Space Telescope Imaging Spectrograph 
(STIS) on board HST showed that the spectrum could be fitted well using a model with a 
mass-loss rate of ~ 10~ 3 M yr -1 . The minimum mass of the system is currently estimated 



to be 120 M ( Hillier et a! 



2HQJ. 



On the basis of photometric and radial velocity variations (in particular the disappearance 
of high-excitation lines like He I, Fe III) and x-ray obs ervations, it appears that n Car is 
a binary system, with an orbital period of 5.52 years (|Damineli. Conti. and Lopes 



Damineli et al. 



200oi 



Ishibashi et al. 



1997 



19991 ) . The mass of the secondary is not known, but 
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it may be less the 30 M (IHillier et adbOQlh . 

Eta Carinae is an extremely enigmatic object on many fronts. However, the morphology 
of its ejected nebula becomes particularly intriguing when we realize that HST observations 
reveal almost identical nebular morphologies in objects of very different mass and evolu- 
tionary history. One of the best known of th ese is the super nova 1987A (SN 1987A) in the 
Large Magellanic Clouds (see, for example, iMcCravl . 120031 for a review). The supernova 



was first observed in February 1987 (and hence 3 years prior to the launch of HST) and was 
immediately classified as a Type II Supernova (SN II; representing the collapse of the iron 
core of a massive star) by virtue of its strong hydrogen lines (coming from the hydrogen-rich 
envelope). For the first time, the detection of neutrino events (formed copiously as matter 
rapidly neutronizes) directly confirmed the association between SNe II and core collapses 
of massive stars. The exploding star itself, SK — 69°202, had actually been observed prior 
to the explosion to be a B3 blue supergiant, with a luminosity of L ~ 1.1 x 10 5 L . Since 
the HST launch, SN 1987A has become a prime target for the telescope, being the nearest 



supernova in modern times. The HST observations have revealed a remarkab 
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e system o f 



Pun 



1997). 



circumstellar rings surrounding the bright center (Fig. 15; iBurrows et al. 
These rings reflect the morphology of material ej ected by the supernov a progenitor a couple 
of tens of thousands of years before the explosion ((Burrows et all Il995i) . This can be inferred 
from the fact that the central ring is expanding at about 10 km s _1 (ICrotts and Heathcote . 
1991 ) and it currently has a radius of about 6.3 x 10 17 cm. While a full explanation for the 



formation of the rings is still lacking, there is very little doubt that what we are observing 
is a bipolar structure, in which the inner ring marks the narrow "waist," while the larger 
rings are somehow "painted" on the bipolar lobes, or possibly mark their edges. The entire 
structure is thus very similar to the one observed in rj Car. 

Bipolar structures have not been restricted only to massive stars, however. Some plane- 
tary nebulae and symbiotic nebulae exhibit morphologies that are almost identical to those of 
r) Car and SN 1987A. Some of the best examples are probably My Cnl8 (Fig. 16; the "hour- 



glass" nebula, usually classified as a pl anetary nebula 



now recognized clS db symbiotic nebula ([Corradi et al. 



, and the "Southern Crab" (Fig. 17), 



2001). 



Planetary nebulae represent the late stages in the lives of stars of about 1-8 M . At that 
phase, the stars eject their outer envelopes exposing the hot cores (T £ 30, 000 K), which, in 
turn, ionize the nebulae causing them to fluoresce. Symbiotic nebulae have at their centers 
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symbiotic binary systems, consisting typically of a red supergiant and a white dwarf that 
provides the ionizing radiation. 

The main question that arises, therefore, is: What is (are) the mechanism(s) that is 
(are) capable of producing such bipolar morphologies in stars of different masses, and dif- 
ferent ages and evolutionary histories? The two main mechanisms that have been proposed 
are: (1) Interacting winds in the presence o f an equatorial to polar density contrast, and 
(2) Magnetic tension of a toroidal field (see iBalick and FrankL 120021 for a more extensive 
discussion) . 

Let us first examin e the interacting winds m odel. The original "interacting winds" model 



for planetary nebulae (|Kahn 



1982; 



Kwok 



19821 ) suggested that old, intermediate-mass stars, 



in the phase of evolution known as the asymptotic giant branch (AGB), first emit a slow 
(~ 20 km s _1 ; of the order of the escape velocity from the AGB star's surface) wind, followed 
by a fast wind (~ 1000 km s _1 ), once the hot and compact nucleus (the AGB star's core) is 



expos ed. The fast wind catches up with the slowly moving material and shocks it. 



Balick 



( 19871 ) proposed that when the interacting winds are allowed to operate in the presence of a 
density contrast between the equator and the pole, a variety of axially symmetric morpholo- 
gies can be obtained. The idea is that, for reasons that will be explored below, the slow wind 
contains a non-spherical density distribution, with material being denser around the equator 
than in the polar direction. Consequently, the (spherically symmetric) fast wind can pene- 
trate more easily at the poles, forming an axisymmetric nebula. Numerical simulations have 
shown that when a range of density contrasts (between the equatorial and polar directions) 



1996 



Frank et al 



1993; 



Icke et al 



1992; 



Mellema 



1995 



be reproduced ( 


Dwarkadas et al. 


Soker and Livio 


1989 


)• 



The second class of models involves the action of a magnetic field. 



The toroidal field component in an outflow from a star is given by (jParksl . Il99lh 

2 



Vr 



rot 



VwiP) 



Rs 
r 



r 



1 sinfl 



(11) 



where Rs, Bs are the stellar radius and the surface magnetic field, V rot is the equatorial 
rotational velocity and Vw{0) is the wind terminal velocity. Consequently, the ratio of the 
toroidal to radial component increases like (r/Rs) at large distances, leading potentially to 
an axisymmetric configuration, by the fact that magnetic stresses can slow down the flow in 
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the equatorial direction (while not interfering with the polar direction). The key physical 
parameters determining the obtained morphology are the stellar rotation rate, the ionizing 
radiation, and the stellar magnetic field. These can be expressed by: Q TOt /Q CTit , F*, and a 



(Chevalier and Luo 



1994; 



Garcia-Segura et al 



1999). Here ^l C nt is the critical (Keplerian) 
angular velocity, f2 crit = (GM S / i?!) 1 / 2 , is the average flux of Lya photons (of order 10 45 - 
10 47 s _1 for typical PNe nuclei), and a is the ratio of the mag netic energy density to the 



kinetic energy density in the wind (e.g., 

B 2 



Begelman and Li 



19921) 



a 



R2 d2 
B S K S 



^pV 2 



\Vw) 



(12) 



■ w MV W W w 
Let us now examine the effects of each one of these parameters. 

For the rotation to have a significant effect and produce a bipolar morphology, the 
star needs to rotate at a significant fraction of its breakup speed (Q TOt /Q cr it ~ 0.5). Un- 
der these conditions, conservation of angular momentum results in a significant focusing of 



the wind towards the equatoria 



Biorkman and Cassinelli 



1993; 



plane, leading to an equatoriall y compressed outflow (e.g., 



Livio, 



1994; 



Qwocki et al 



1994). The concomitant equator 



to pole density contrast leads to bipolar morphologies. 

Ionization does not, in itself, produce a bipolar morphology. Rather, ionization fronts tend 
to excite ins t abiliti es (similar to the Rayleigh- Taylor instability or the instability discussed 



by 



Vishniac 



1983), which, in turn, produce finger-shaped structures and dense "knots." 



Observations of a number of relatively nearby planetary nebulae, and in particular of the 



"Helix" nebula, re veal that such dense knots are probably very common ([O'Deil et al 



Speck et al 



2002). 



2003; 



The toroidal magnetic field is carried by the fast wind and it can, in principle, pro- 
duce a bipolar morphology even if the slow wind is spherically symmetric. Basically, as 
the magnetic energy density of the shocked wind becomes larger than the thermal energy 
density, the flow becomes (due to the increasing importance of the toroidal component) 
bipolar, with an increased c ollimation as the value of a is increased. Numerical simulations 
( Garcia-Segura et a/.Ul999j ) show that qualitatively, the morphologies obtained for different 



combinations of n rot /f2 cr i t and a are as shown in Fi g. 18. The minimum fie 



produce magnetic shaping was found to be of order ( Chevalier and Luo 



Dmm 
B S 
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M 
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FW 
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199J) 



0.1 



(13) 
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where Vpw is the velocity of the fast wind. 

From eq. 13 we see that if the star rotates too slowly, the minimum required field may 
be unattainable. Similarly, if the star does not rotate at a significant fraction of its breakup 
speed, an equatorially compressed outflow is not formed. The question of: What is the 
mechanism that produces highly bipolar outflows? can therefore be reduced to: What 
causes the star to rotate close to breakup? or: What can produce a strong density contrast 
between the equatorial and polar directions in the slow wind? 

One obvious possibility is: binary companions! 

Companions to the central star can act in several ways to aid in the formation of bipo- 
lar morphologies: (1) For binaries that were initially relatively close (separation less than 
~ 1000 R©), so that the primary could fill its Roche lobe (the critical potential surface be- 
yond which mass transfer onto the companion occurs) during the asymptotic giant branch 
(AGB) phase, an unstable mass transfer ensues. As a result, the companio n and the AGB 



star's core start spiralling-in inside a common envelope (see, e.g., a review by llben and Livio . 



19931 ). This has two effects. First, the envelope of the primary can be spun-up to angular 



velocities of the order of 



Oat n J M c \ (M env \ 1 (K 2 g \ 1 / a \V2 



n CIit - 1 [omm g )[m q J [ol) [rs) ' (14) 

where Mc is the companion's mass, M env is the mass of the giant's envelope, M$K 2 R 2 S 
is the star's moment of inertia and a is the initial separation between the giant and the 
companion. Equation 14 shows that even brown dwarf (sub-stellar) companions can bring 
the envelope close to critical rotation. Second, since the envelope mass is ejected (due to 
orbital energy deposition) primarily close to the orbital plane (because angular momentum 
is also deposited into the envelope), the common envelope phase can generate an equator-to- 
pole density contrast. Hydrodynamic simulations of common envelope evolution reveal that 
durin g the late stages a bout 80% of the mass is ejected within 30° o f the b inary orbital plane 



e.g. 



Rasio and Livio . 



1996 



Terman. Taam. and Hernquist 



1994 . 



19951 ). One can expect 



that due to cooling, the mass will sink even more toward the orbital plane at later times. 

Another possibility for the presence of a higher equatorial density, in principle at least, 
that does not involve binary companions, is the inner rim of the protostellar disk. If the 
outer part of the protostellar disk survives till late stages in the stellar life (which may not 
be difficult in the case of massive, short-lived stars), then the fast wind could interact with 
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the inner rim of this disk (jPringle, 198 91) . In a few planetary nebulae (for exam ple, the "Red 



Rectangle" and the "Egg Nebula," iBond et al 



1996 



Thompson et al 



19971 respectively), 



HST observations reveal the presence of relatively large disks, similar to the ones observed 
in young stellar objects. 

In at least some symbiotic nebulae (e.g., M2— 9), the bipolar morphology may reflect the 
action of the white dwarf companion. The white dwarf accretes from the wind of the AGB 
star, an accretion disk is formed, and t he disk powers a mildly collimated fast wind, which i n 
turn produces the bipolar morphology ([Livio and Sokerl . l200Tt ISoker and RappaportL 120011 ). 

The conclusion from this discussion is that several mechanisms are capable, in principle 
at least, to produce the observed bipolar morphologies. Different mechanisms may be op- 
erating in different systems. In some cases, we can look forward to the future and expect 
more definitive answers to emerge. For example, in SN 1987A, the supernova blast wave will 
eventually hit the entire inner ring, and the luminosity tha t will be generated by this inter- 
action will illuminate the entire SN vicinity (jMcCravl . 120031 ) . A few brightening spots, where 
the blas t wave ha s alrea dy hit protrusions on the ring, have been observed by HST (Fig. 19; 



see also 



Panagia, 



20021 1 . Generally, a transmitted shock at normal incidence is expected 



to propagate into the ring with a speed of V ring ~ (w /^ing) I'biast (where n , n ving are 
the number densities of the circumstellar matter and the ring, respectively, and Vbi as t is the 



blast wave velocity). For SN 1987A, V Uast ~ 4000 km s" 1 



, n D ~ 150 cuf 



, f t-rmg 



10 4 cm 



giving V r - mg ~ 500 km s 1 . Eventually, the Ha flux from the ring is expected to be more than 
30 times higher than today (Fn„ ~3 x 10~ 12 erg cm" 2 s -1 ), and even brighter in UV lines 



fhuo. McGrav. and Slavin. 



19941 ). This ionizing flux will turn the circumstellar matter into 



an emission nebula, thus revealing its distribution and velocity field, and hopefully allowing 
for a reconstruction of the mass-loss history of the system. In all of this, HST will provide 
a front seat view. 

The collapses of very massive stars are believed to produce stellar-mass black holes. In 
some cases, these collapses produce the most spectacular explosions in the universe since the 
"big bang" — Gamma Ray Bursts. The Hubble Space Telescope has played an important role 
in the attempts to understand the nature of these dramatic explosions. However, much more 
massive black holes are also produced in the universe quite commonly, and the questions of 
how they form and evolve have intrigued astronomers for decades. The initial answers had 
to come from large collections of stars like galaxies and stellar clusters. 
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Given the unprecedented resolution that HST provides, very crowded fields, like the 
centers of galaxies, were obvious targets. In particul ar, the early suggestive evidence that 
galactic centers harbor supermassive black holes (see iKormendv and Richstond . Il995l for a 
review), virtually ensured that centers of galaxies will be observed extensively with HST. 
Indeed, these observations were carried out and proved to be extremely fruitful. 



V. BLACK HOLES FROM SUPERMASSIVE TO STELLAR 



The idea that active galactic nuclei (AGNs) are powered by accretion at high rates (up to 
a few solar masses per year) onto supermassive (10 6 -10 8 M m) central black holes has been 
around for a long time (e.g., IZeldovich and Novikovl . Il964l ). Furthermore, it has long been 



realized that collapse into a black hole may be inevitable in the center of man y dense stellar 
environments, on ce the central potential gets sufficiently deep (for example, Lvnden-Bell 



1969; 



RCC! 



19781 ) . Consequently, one might expect supermassive black holes to reside at the 



center of most galaxies. 

In addition, searches for the host galaxies of quasars (thought to represent accret- 



ing supermassive black h o 



neer ing works of 



and 



Kristiarj 



es) h a ve been going on fo r a long time, following the pio - 
19731 ) Iwvckoff et all l|l98(j h Irlutchings and Camobelll (|1983l ). 



Hutchings. Janson and Nefa ()1989l ). The ground-based studies not only detected host 



galaxies, but also sugges t ed that some of the hos ts have experienced tidal interactions 



(jBoroson. Oke and Green . 



1982 : 



Kukula et al. 



1996). 



These tentative detections and suggestions have turned into re ality with the spec 



tacular, high-resolution images of quasars obtained wi t 



Bahcall. Kirhakos. and Schneider 



1994 : 



Kirhakos et al 



hHST (e.g., 



Bahcall et al. 



1997 



19991 ) . The HST images revealed 



unambiguously that nine radio-loud quasars reside either in bright elliptical hosts or in 
strongly interacting systems. These findings probably indicate that interactions with neigh- 
boring galax ies play an im portant role in the fueling of the central black hole in active 
galaxies (see 



Krolik 



1999ah . 
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A. Search techniques for supermassive black holes 



The actual searches for s upermassive black holes rely pr imarily on stellar- and gas- 



dynamical evidence (see, e.g., 



Kormendv and Richstone . 



1995 



for a review). In particular, 



the idea is to unambiguously show that the mass-to-light ratio, M/L, increases toward the 
galactic center to values that are difficult to accomodate with other types of stellar popu- 
lations. Ideally, one would want to follow this up with the detection of relativistic speeds, 
but at present even HST cannot resolve orbits at a few gravitational (Schwarzschild) radii 
(Rs = 2GMbh/c 2 ; where Mbh is the black hole mass, G is the gravitational constant and c 
is the speed of light). 

Generally, the basic principle behind the early stellar-dyamical search techniques can 
be explained by the following, simplified picture. Taking the first velocity moment of the 
collisionless Boltzmann equation gives 



M(r) 



G G 



dlnpt dlnai 



dlnr dlnr 



On 



cr: 



at 



(15) 



where M(r) is the mass enclosed within radius r, V TOt is the rotational velocity, a r , <j e , 
0$ are the components of the velocity dispersion, and pt is the density of the tracer stars 
being observed (usually assumed to be proportional to the volume brightness). A direct 
measurement of V TOt and a, therefore can determine the central mass. To actually use eq. 15, 
however, the ranges of unprojected quantities (e.g., V rot , a) need to be derived, and various 



techniques to achieve that have 



1988; 



Gerhar 



d 



1993; 



Kormendv 



jeen developed (see for ex ample, iDressler and Richstone , 



1988; 



van der Mare l[l99j). 



The more recent search techniques fit axisymmetric, three-integral dynamical models of 
the galaxy (using the line-of-sight velocity distribution), to the observed light distribution. 
Basically, the orbits in the (R, z) plane in realistic galactic potentials are often found to 
have in addition to the two integrals of motion E (the energy) and L z (the ^-component 
of the angular momentum, where z is the symmetry axis), a third integral, J3, that can be 
associated with the approximately conserved total angular momentum, L. 

The dynamical model assumes axisymmetry and an inclination for the galaxy and first 
determines the optimal density distributrion that is consistent with the surface brightness 
distribution (making certain assumptions about M/Ly). Then, a central point mass is 
added, and the potential calculated. Orbits that sample the phase space of the three integrals 
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of motion are then calculated, and the data of the full line-of-sight velocity distribution are 
fitted to the models. 

Since galactic gas is affected also by forces which are non-gravitational (e.g., radiation 
pressure), stellar kinematics are considered more secure than gas dynamics in black hole 
searches. Nevertheless, in a few cases, the presence of the black hole may be revealed by 
gas-dynamical searches. The prototype of this technique is provided by the radio galaxy 



M87. Early stellar-dynamical observations revea 



rise inward, to r ~ 1/(5 (e.g., 



Lauer et al. 



1992: 



ed that the ve l ocity dispersion continues to 



Sargent et a/l Il978[ ). However, due to the 



expected anisotro py in the velocity dispersion , models without a black 
obser vations (e.g., 



Binnev and Mamon . 



1982; 



Dressier and Richstone 



role could also fit the 

1 

199(1' 



van der Marell . 



19941 ) . In particular, the last author found a ~ 400 km s 1 at r <,0"5, which could never- 



theless be fitted with anisotropic models that do not include a black hole. 

A more definitive answer, however, came in this case from the gas. High-resolution 
HST observations of the nucleus revealed the presence of a small gas disk (about 20 pc in 
radius), with a major axis perpendicular to the optical jet (Fig. 20). Spectra taken at two 



emission lines separated by 2V = 916 km s 1 f)Ford et al 



(Ford et al. 


1994: 


Harms et al.. 


1994) 



an inclination angle of the disk of ~ 42° (implied by the observed axis ratio), the observed 
velocity amplitude (if interpreted as a circular Keplerian motion) corresponds to a dark 
m ass of Mrh ~3x 1 9 at the center of M87. These results have been further confirmed 



by 



Macchetto et al 



( 19971 ). making the black hole in M87 the most massive observed so far 
(with a relatively secure mass determination). I should caution that an airtight case for 
circular Keplerian motion is still to be made for M87. Incidentally, the best case for a black 
hole based on gas dynamics is for the modest active galactic nucleus NGC 4258. In that 
case, Very Long Baseline Array (VLBA) observations of high- velocity masers show a rotation 
curve that can be fitted remarkably well with V rot (r) ~ (832 ± 2)(r/0.25 pc)" 1 / 2 km s _1 , 
implying a black hole mass (mass interior to 0.18 pc) of Mbh — 4.1 x 10 7 M Q . 



B. The Mbh a relationship 



Using three-integral models and HST spectroscopy, 



Gebhardt et al. 



(2000J) were able to 



determine black hole masses in 16 galaxies. Adding to these two galaxies with maser mass 
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determinations (NGC 4258 and NGC 106 8), six galaxie s with black hole masses d etermined 



by gas dynamics, and our ow n galaxy ([Genzel et al 



( Dressier and Richstone . 



2000; 



Ghez et al 



19981) and M31 



1988), they were able to discover a tight correlation between the 



black hole mass and the velocity dispersion in the galactic bulge. Specifically, they found 
that 



Mi 



BH 



1.2 x 10 8 M P 



o 



3.75 



200 kms"V ^ 
where a e is the line of sight aperture dispersion within the half-light radius R e . The ob- 

tained relation is much tighter (Fig. 21, right) than a previously determined relation (e.g., 



Magorrian et al. 



1998) between black h ole mass and the bulge lumi nosity (Fig. 21, left). 



Based on a smaller sample (12 galaxies) . iFerrarese and Merrittl (j2000l ) found independently 
a Mbh oc a a relation, with a somewhat steeper slope of a ~ 4.8. 

Clearly, the observed tight M B h — cr correlation strongly suggests that the formation and 
evolution of the bulge and of the black hole are causally connected. The precise nature of 
this connection, however, is still a matter of consider able uncertainty. A relatively simp 



theoretical model for the relation has been suggested bv lAdams. Graff and Richstonel ()200lh . 



In this model, a slowly rotating isothermal sphere (with a seed central black hole) collapses 
to form the bulge (the dark matter is assumed to move in tandem with the baryons). The 
density distribution is assumed to be of the form 



p(r) 



(17) 



2nGr 2 ' 

where c s is the speed of sound. Note that for dissipationless collapse the velocity dispersion 
is roughly given by o 2 ~ 2c 2 s . The region is assumed to rotate rigidly (due to effective 
tidal torques) at an angular speed fl, and the specific angular momentum is j = r^Qsm 2 9, 
where is the initial radius and 6 is the polar angle. For zero-energy orbits, the pericenter 
distance in the equatorial plane is therefore 



V 



r 



(GM) 3 Q 2 



(18) 



2GM 2GM 32c 8 s 
where in the last equality we used the fact that for the assumed density distribution M{r) = 
2c 2 s r/G. For material to be captured by the black hole we need p < 4i?s, where Rs is the 
gravitational (Schwarzschild) radius Rs = 2GM/c 2 . Using eq. 18, the condition p = 4R S 
therefore reads (assuming that the capture condition defines the black hole mass) 

16cf 4 , 



Mi 



BH 



GcVt 
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Gctt 



a 



(19) 



in good agreement with the observed relation (eq. 16). 

Other considerations result in similar expressions. For example, in a protogalaxy modeled 
as an isothermal sphere of cold dark matter, with p(r) = a 2 /27rGV 2 , with a fraction / gas in the 
form of gas, a central, accreting black hole will generate an intense wind outflow. The black 
hole itself may be assumed to form by coherent collapse before most of the bulge gas turns 
into stars. If the black hole radiates at the Eddington luminosity, Ledd = 47tcGMbh/ k (at 
which gravity is balanced by radiation pressure; where k is the electron scattering opacity), 
and a fraction / out is deposited into kinetic energy of the outflow, then a shell of swept-up 
material will be moving outward at a speed 



v; 



out 



87T 2 G/ out LEDD 



1/3 



(20) 



The condition that the shell would escape, and therefore, that the black hole would unbind 
the bulge gas, requi res V out > a. This implies that the black hole mass is limited by 



(Silk and Ree, 1 



19981 ) 



Mi 



BH 



similar to the relation found by 



,w/ ou t r?v 



o" 



(21) 



Ferrarese and Merritd (|2000j ) 



Semi-analytical, hierarchical galaxy formation models (see Section VI) , in which galaxies 
form by merging halos (and merging central black holes), with simple prescriptions for gas 
cooling, star formation, and feed back from supernovae, also ten d to produce scalings of 
the form Mbh ~ c 4 (for example. lHaehnelt and Kauffmannl . I2000T ) . Broadly speaking, this 



relation can be traced to the facts that: (i) In mergers, the black hole mass scales with the 
halo mass, (ii) o ~ p 1//6 Mj|^, and (iii) M ha i scales like p~ 2 in typical cold-dark- matter 
cosmologies (really like p~ 2 /( 3+ri ), where n ~ —2 is the slope of the dark matter fluctuations 
spectrum). Combining (i)-(iii) gives M B h ~ cr 4 - 

Somewhat more exotic scenarios, which involve the accretion of collisional dark m atter 
(invoked to make galactic halos less dense; for example. ISper gel and SteinhardtL 120001 1 . also 
produce black hole masses which scale roughly with a 4,5 ( Ostriker . I2OO 



C. Intermediate- mass black holes? 



An interesting question is related to the range of black hole masses over which the ob- 
served relation (16) applies. Most recently, high-spatial-resolution spectroscopy with HST, 
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(Gerssen et al 



2002, 



2003 



van der Marel et al. 



20021 ). of the central part of the globular 



cluster M15, r evealed evidence 
mass (see also iBaumgardt et al. 



'or the existence o 



2002 



Dull et al. 



a de nse, central concentration of dark 



2002 ). The interp reta tion of the nature 



of this dark mass has been somewhat controversial. 



Dull et al. 



(2002) and 



Baumgardt et al 



(2002J) have shown that the sharp rise in M/L that was observed toward the center of M15 
could be explai ned by a centra l conc entration of neutron stars and massive white dwarfs. On 
the other hand, iGerssen et al\ (J2003J) argued that if one were to allow for the fact that a large 
fraction of the neutron stars born within the cluster cannot be retained (due to "kick" veloc- 
ities of a few hundred km s _1 at birth), the observed M/L might require the existence of a 
central black hole of mass Mbh = l- 7l? 'I x 10 3 Mm. Interestingly, the derived ma ss fits quite 



well onto the M B h — a relationship (|Gerssen et al 



2002; 



Tremaine et al. 



2002J). Similarly, 



evidence for a central black hole with a mass of ~ 2 x 10 4 M^has 



cluster Gl in the Andromeda galaxy M31 (|Gebhardt. Rich, and Ho 



Deen found in the stellar 



200 21). The latter blac k 



hole also fits the Mbh — cr relation. In an entirely independent work (jD'Amico et all 12002). 
pulse timing observations determined positions for five millisecond pulsars in the cluster 
NGC 6752 with a 20 mas accuracy. Three were found to have line-of-sight accelerations 
larger than the maximum value predicted on the basis of the central mass density derived 
from optical observations. The measured accelerations thus provided dynamical evidence for 
a central mass-to-light ratio of M/L £ 10. All of these findings, if confirmed, are extremely 
exciting, since the implied masses make these black holes intermediate between the stellar 
black holes (with Mbh ~ 10 M ) and the suppermassive ones (Mbh ~ 10 6 -10 9 M ). 

The potential existence of a class of intermediate-mass bl ack holes has also been suggeste d 
following the discovery of ultraluminous x-ray sources (e.g.. IColbert and Mushotzkvl .ll999). 
The latter are x-ray sources outside the nuclei of external galaxies, with luminositi es in excess 
of 10 39 erg s _1 . These sources were originally discovered by the Einstein satellite (jFabbiano, 



1! 

ries 



b ut have been found in large numbers by the ROSAT and Chandra observato- 



e.g.- 



Colbert and Mushotzkv. 



1999; 



.Tel tern a et al 



2002; 



Lira. Johnson, and Lawrence , 



2002T ). More recently, however, it has been pointed out that although the ultraluminous 
x-ray sources may form a heterogeneous class of systems, the most likely explanation for the 
majority of them is that they constitute the high-luminosity tail of the ste l lar-mass black- 
hole binary distribu tion (e.g., 



200lh . In particular, 



King et a 



Podsiadlowski et al 



2001 



3) 



Podsiadlowski et al 



2002; 



Roberts et al. 



(2002J) have shown that in binary systems in which 
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the donor star becomes a giant and the evolution is driven by the nuclear evolution of the 
hydrogen-burning shell, luminosities that are potentially as high as ~ 10 41 erg s -1 can be 
obtained. The existence of ultraluminous x-ray sources by itself, therefore, may not imply 
the existence of intermediate-mass black holes. 

The potential existence of intermediate -mass black holes also led to a new scenario for 
the formation of the supermassive ones ((Ebisuzaki et all 120011 ). The idea is that first, 
intermediate-mass black holes form in young clusters, due to runaway mergers of massive 
stars. The massive sta rs are assumed to sink to the cluster center due to dynamical friction, 



on a timescale of (e.g., 
1.17 r 2 a 



Binnev and Tremaine . 



mi 



log A Gm 



2.7 x 10 7 



n 



-1/2 



M 



1/2 



m 



1 pc J \ 10 pc 



10 6 Mr, 



20 M m 



yr 



(22) 



where log A is the Coulomb logarithm, a is the velocity dispersion, r is the distance to the 
cluster center, is the half-mass radius and M and m are the masses of the cluster and 
the star, respectively. During the same time that the intermediate-mass black holes are 
forming, the host clusters themselves sink to the galactic center (again due to dynamical 
friction), evaporate, and deposit their black holes. The latter can then form black hole 
binaries which merge due to gravitational radiation, eventually leading to the formation of 
supermassive black holes. I should note that this is only one path out of the many that ha ve 
been suggested over the years for the formation of supermassive black holes (jReesl . Il984l ). 

The installation of HST's Advanced Camera for Surveys (ACS) in March 2002, with 
its superior sensitivity (by a factor of about 3-5 over the previous Wide Field Planetary 
Camera 2) and increased resolution (in the high-resolution channel), combined with X-ray 
observations by the Chandra and XMM-Newton Observatories, and infrared observations 
of dust-obscured active galactic nuclei by the Space Infrared Telescope Facility (SIRTF) 
promise that the study of black holes in clusters (if they exist) and in galactic centers is only 
beginning. In particular, future observations will help clarify the relation between galaxy 
evolution in general, and the formation and evolution of the cental black holes. 

Observations with HST have provided important insights not only in the study of super- 
massive black holes, but also in researches related to the probable formation of stellar-mass 
black holes — the study of gamma-ray bursters. 
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D. Gamma-ray bursts 



Gamma-ray bursts were first detected by the military Vela satellites in 1967. These are 
short flashes of gamma-rays (typically with a peak energy around 100 keV), lasting between 
a few milliseconds and tens of minutes. Already in the early 1990s, the BATSE experiment 
on board the Compton Gamma- Ray Observatory demonstrated that the rate of gamma-ray 
bursts is about 1-2 per day, and that they are distributed isotropically in the sky. The 
data from BATSE also showed that gamma-ray bursts (GRBs) come in two distinct classes: 
(i) short-duration (£, 2 sec) with hard sp ectra, and (ii) long duration with softer spectra 



(Fishman, 



2001 



Kouveliotou et al. 



199c 



Due to the relatively poor localization capabilities of gamma-ray detectors (which did 
not allow identification of the sources in other wavebands), two thousand bursts have been 
detected before it became possible to determine the distance to the bursts. The isotropic 
sky distribution argued for either a cosmological origin or an extremely local one (e.g., the 
Galactic halo). This situation changed dramatically since the operation of the Italian-Dutch 
BeppoSAX satellite in 1997. BeppoSAX was able to determine the positions of bursts to 
within arcmi nutes. The d i scovery of rapidly declining "afterglows " in the X-ray, optical anc 



radio bands (Costa et al 



1997 



Frail et al. 



1999: 



Metzger et al. 



1997: 



van Paradiis et al. 



19971 ) allowed for redshift determinations that immediately placed GRBs at cosmoligical 
distances. I should note that to date, only afterglows of the long-duration bursts have been 
observed. Thus, it is not even clear if the short bursts produce afterglows (the best limit to 
date, for the burst GRB 0205 31, failed to detect an afterglow candidate [down to V ~ 25] 



about 20 hours after the burst; 



Salamanca et al 



2002J). Typically, the afterglows decay with 



time as t~ a , with a in the ranges 1-2. This behavior was predicted by a mode l in w hich a 
fireball is expanding into a homogeneous external medium ([Meszaros and Reed . 119971 ) . 

At the time of this writing redshifts have been determined to more than two dozen GRBs, 
and they usually lie in the z = 0.5-1.5 range (although redshifts as high as 4.5 have been 
recorded). The observed fluxes imply energies of up to 10 54 ergs for isotropic emission. 
However, there is increasing evidence, in the form of kinks in the afterglow light curve, and 
in polarization detected in a few bursts, that gamma-ray bursts are in fact collimated into 
narrow jets. If the observer's line of sight is within the jet solid angle, Qj, then as long as the 

— 1/2 

Lorentz factor 7 satisfies 7 ^ fL , the light-cone is within the jet boundary. However, as 
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the jet decelerates, 7 eventually drops below fl • . Consequently, the (transverse) emitting 
area starts to grow more slowly (as r^nj 1 ^ 2 instead of (r \\ H) 2 ), resulting in a break (faster de - 
cay) in the light curve (|Livio and Waxman . 



2000; 



consistent with observations in GRB 990123 rtFruchter et al 



Meszar os and Rees l 



1999; 



1999; 



Rhoads, 



Kulkarni et al. 



1997), 



1999). 



When collimation is taken into account, the average total energy of GRBs is e stimated to 



be of the order of 2 x 10 51 ergs ([Frail et al 



2001 



Panaitescu and Kumar 



20011). 



Based on the energetics, and the fact that the production of astrophysical jets typically 
relies o n the collim ation and acceleration provided by an accretion disk around a compact 



object (jLivio . 



2000), the most popular models for GRBs involve the formation of stellar- mass 



black h oles, surrounded by a debris torus from which mass accretion onto the central object 



occurs ( Meszaros , 



2002). The two most likely progenitors to produce such a configuration are 



the core collapse of massive stars ([Paczvn ski 



1998 



Wooslev . 



1993), or the ongoing 



of neutron star-neutron star or b lack hole- neutron star binaries (e.g., 



Goodman . 



1986 



Paczvnski 



merger 



Eichler et al 



1989; 



1986). In both cases the main source of energy is gravitational, 



even though tapping into the 



(Bla ndford and Zna jek. 



12Z7J 



arge) spin ene rgy of the black hole is also possible in principle 



Krolik . 



1999bl ). A schematic describing how a GRB and its 



afterglow may be generated by internal and external shocks (respectively) in a collimated 
jet arising from a stellar collapse, is shown in Figure 22. 

There are several pieces of evidence suggesting that the long duration GRBs may indeed 
be associated with stellar collapses. First, in merging neutron stars (or a black hole and 
a neutron star) the associated timescales (dictated at that point by the emission of grav- 
itational wave radiation) are probably too short to produce multi-second-long bursts. On 
the other hand, core collapses that lead to a black hole and an accretion disk are naturally 
associated with longer timescales. Second, there exists quite strong evidence that at least 
some GRBs are associated with supernova explosions. 

In particular, the supernova SN 1998bw (of the relatively rare type Ic, caused by the 
explosion of a massive star th at had rid itself of its hydrogen envelope prior to exploding; 



Filippenko and Sargent 



1986) was found to be approximately coinci dent in both positio n 



and time with the relatively weak GRB 980425, at redshift z = 0.0085 (jGalama et a/.L Il998). 
Furthermo re, the supernova ligh t curve was found to be consistent with the formation of a 



black hole (Iwamoto et al 



1998). In a few other bursts (e.g., GRB 011121) the optical and 



near infrared afterglow light curve exhibits a "bump" with a time delay and amplitude that 
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are consist ent with resulting fro m a supernova explosion occuring simultaneously with the 
GRB (e.g. 



Greiner et al. 



2001). 



The contributions of HST to this exciting field have been in a few areas. 
First, by resolving the host galaxies and being able to pinpoint the GRB on the host, 
HST has shown unambiguously that at least some GRBs are not associat ed with galactic 



nucle i . In fact, they som etimes occur far from the nucleus or in spiral arms (lAndersen et al 



2002; 



Sahu et 



Fruchter et al 



al 



1997). Second, the Hubble Space Telescope has shown ((Fruchtei, 



2002; 



1999T ) that the colors of the hosts of GRBs corresond to the bluest colors 



observed for galaxies in the Hubble Deep Fields (see Section VI). This means that GRBs 
with afterglows occur preferentially in galaxies with high star formation rates — a finding 
that is consistent with GRBs being associated with collapses of massive stars. Furthermore, 
the fact that in almost all cases the GRB's position was found to be within the extent of 
the rest-frame ultraviolet image (where star formation is intense) of the hosts, also argues 
in favor of the "collapsar" model for the long-duration bursts. Neutron star binaries, on the 
other hand, are born with "kicks" of as much as a few hundred km/sec, which could drive 
them out of the star-forming regions. 

The Hubble Space Telescope also helped to firm the association of at least some GRBs 
with supernova explosions. In GRB 011121, for example, HST detected an intermediate- 
time flux excess ("red bump") that was redder in color relative to the GRB afterglow. 
This "bump" could be well described by a r edshifted Type Ic supernova that explode d 



approximately at the same time as the GRB (Bl oom et al 



2002; 



Garnavich et al 



2003). 



Near-infrared and radio observations of the afterglow further provided evidence for extensive 
mass loss (M ~ 2 x 10 -7 M yr _1 ) from the massive stellar progenitor of GRB 011121 



(Pri ce et al 



2002). 



To conclude this topic, it is very likely that every time a GRB goes off, a black hole is born. 
Gamma Ray Bursts therefore offer, in principle at least, a direct measure of the massive-star 
formation rate in the early universe. Future observations with the HETE-2 spacecraft and 
with the Swift multi-wavelength GRB afterglow mission (equipped with 7-ray, X-ray and 
optical detectors), as well as follow-ups with HST and other observatories from the ground 
and space, will hopefully reveal the true nature of these most dramatic explosions (including 
the short- duration bursts). 
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VI. DEEP FIELDS AND A BRIEF COSMIC HISTORY 



One of the major goals of observational cosmology is to understand the processes involved 
in galaxy formation and evolution. A full theoretical treatment of this problem requires an 
understanding of the primordial density fluctuations, of the formation of dark matter halos, 
of the dissipation that occurs in the cooling gas, of the processes involved in star formation, 
of the feedback between stellar explosions and the interstellar medium, of galaxy mergers, 
and of the interactions of galaxies with intergalactic gas. 

Most current models for the formation of structure in the universe assume that dark 
matter halos build up hierarchically, with the assembly being controled by the cosmological 
parameters, the power spectrum of the density fluctuations, and the nature of the dark 
matter itself. The build-up of the stellar mass within galaxies is a secondary process, in 
which gaseous dissipation, the more intricate physics of star formation, and feedback, all 
play a role. 

Since galaxy formation involves non-linear processes on many scales, the problem has 
been addressed th e oretic ally mainly using large-scale iV-body hydrodynamic simulations 



(.e.g. 



Pearce et al. 



20011 1. On the observational side, HST, the Chandra Observatory and 



the Keck and VLT telescopes have allowed for an unprecedented view of the high-redshift 
universe, probing galaxies from infancy to old age. 

Short of the awsome pillars of dust and molecular gas revealed by the HST images of the 
"Eagle Nebula," the best known images that Hubble has produced are those of the "Hubble 
Deep Fields" (HDFs). In fact, the first HDF and the follow-up observations of the same 
field by other observatories probably represent the most concentrated research effort ever in 
astronomy into what was previously a blank piece of the sky! 

Observations taken shortly after the first servicing mission, which restored the HST optics 
capabilities, demonstrated that HST could resolve distant galaxies spectac ularly well. In 



parti cular, observations of the relatively distant c lusters CL 0939 +4713 (jDressler et al 



19941 ) and the cluster around the radio galaxy 3C324 (jDickinso n. 1995) revealed faint galaxies 



of small angular sizes, the morphologies of which were entirely inaccessible from the ground. 

The idea of the Hubble Deep Field North — ten days of one continuous observation of a 
field 2.6 arcminutes on the side in December 1995 — was conceived by the then Director of 
the Space Telescope Science Institute, Bob Williams. Williams decided to use his Director's 
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Discretionary Time on HST to produce the deepest image of the universe in optical/UV 
wavelengths. 

The northern field itself (Fig. 23) was carefully selected so as to be in HST's northern 
continuous viewing zone (CVZ; without interference by Earth occultations), to be free of 
bright stars, nearby galaxies and radio sources, and to have relatively low Galactic extinction. 
The southern field (Fig. 24) was selected so as to include a quasi-stellar object (QSO) that 
could be used to study absorption systems along the line of sight. The HDF-N observations 
were taken in December 1995 and the HDF-S in October 1998. 

The filter selection for the observations was driven partly by the desire for depth and 
color information (to identify high-redshift galaxies by their Lyman-break) , and partly by 
practical considerations involving scattered light within the telescope. Accordingly, images 
were taken in four broad-band passes, spanning a wavelength range from about 2500 A to 
9000 A. 

In order to achieve a higher resolution than the detectors' pixel sizes, the pointing of the 
telescope was dithered (shifted slightly) during the observation, thus ensuring that images 
were recorded on different pixels. A "drizz ling" (subpixel linear reconstruction) technique 



was developed by 



Fruchter and Hook! (J2002J). Details on the observati onal techn i ques a nd a 

( 19981 ) and 



review of some of the results can b e foun d in the excellent reviews by 



Ferguson 



Ferguson. Dickinson, and Williams! ((2000). 



As soon as the HDF-N image was obtained, it was obvious that a new era in astronomy 
has begun. The image revealed some 3000 galaxi es of different colors, sh a pes, a nd sizes. 



A broad summary of key results can be found in iLivio. Fall, and MadaiJ (|l998f ) . Here I 
will concentrate only on two main topics, to which HST's contribution has been crucial: 
(i) galaxy sizes and morphologies (and their implications for galaxy formation and evolution), 
and (ii) the global, cosmic star-formation history. The dramatic discovery of a particular 
supernova in the HDF-N will be discussed in Section VIII. 



A. The morphology and sizes of high-redshift galaxies 

One of the key questions in galaxy formation and evolution is the importance of interac- 
tions and mergers. In particular, cold-dark-matter and dark-energy dominated models have 
as one of their natural consequences a higher merger rate in the past, since in these models, 
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toda y's Hubble Sequence galaxies have been assembled via a process of hierarchical m ergers 



[e.g., 



Cole et al. 



2000; 



Kauffmann et al 



1999; 



Somerville. Primak. and Faber . 



2001). 



In the local universe, the luminosity density is dominated by ellipticals, lenticulars and 
spirals, even though dwarf ellipticals and dwarf irregulars (all of low luminosity) dominate 
by number. The HDF-N allowed for the fi rst time for a morpho logical classification of 
galaxies down to a magnitude of / = 25 (e.g.. lAbraham et q/.l ll996). These studies showed 



that the fraction of irregular, multiple- component, and peculiar- looking galaxies was indeed 
considerably higher (~ 40%) than expected from a direct extrapolation of the numbers 



at z ~ 0. The trend of a rising fraction of 



already noted in HST's Medium Deep Survey ([Griffiths et al 



aint systems with irregular morphology wa s 



1994 



Windhorst et al. 



1995), 



although there, the survey reached only down to / ~ 23. 

One worry one might have when examining the morphologies of galaxies at high redshift 
is that optical images reflect, in fact, the UV rest frame images of the galaxies. Since 
the UV light traces, in particular, pockets of intense star formation, one m ight expect a 



more irregular appearance for high- redshift galaxies (IGiavalisco et all Il99fif) . A study of 
the redshift distribution of galaxies with 17 < / < 21.5 (jlm et all Il999l ). showed indeed 



the irregular/peculiar class to consist of both low- redshift dwarfs and intermediate redshift 
(0.4 £j z £j 1) galaxies that are by themselves a mix of starbursts (of relatively low mass) and 
more massive, i nteracting galaxies. 



Nevertheless, 



Brinchmann et al. 



( 19981 ) and 



Abraham et al. 



( 1999a| ) have shown that the 



redshifting of UV wavelengths into the optical is not sufficient to explain the preponder- 
ance of irregular morphologies. Further confirmation that the rapid rise with redshift in 
the fraction of galaxies with irregular /peculiar morphologies is real, came from the Near 



Infrared Camera and Multi O 



observations ( Dickinson 



pie ct Spectrograph (NICM OS) observations. The NICMOS 



2000s 



Thompson et all 11999^ ) have demonstrated that with very 



few (interesting by themselves) exceptions, the irregularities seen in the optical (WFPC2) 
image persist in the infrared image (Fig. 25), even at redshifts as high as z ~ 3. Since 
the NICMOS observations reflect the rest-frame optical light (up to z ~ 3), the absence 
of any significant changes in the morphology (in the majority of galaxies) proves that the 
peculiarities are not artifacts of wavelength shifting. 

A morphological trend that may be indicated by the HDFs is a significant de cline in the 
fraction of barred spirals with redshift (for z^,0.5; e.g. 



Abra ham et al. 



1999b). However, 
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since bars are more difficult to detect at bluer rest wavelength s , optical surveys of galaxies 



at hi gh- z may be biased against finding bars ((Eskridge et al. 



2000; 



van den Bergh et al. 



20021 ). The precise reason for this trend (if confirmed) is not known, but it may be related 



either to the state of the disk (the development of a bar instability requires a cold disk), or 
to the timescal e needed to form long-lived bars (perhaps through episodic growth, aided by 



spiral patterns; 



Sclwoo 



d. 



2000). 



Another extremely interesting result to have come out of the HDFs is related to 
the sizes of galaxies. When examining the HDFs, one is immediately struck by the 



small angular diameters o 



Bouwens. Cavon. and Silk 



the faint galaxies. In particular, No- Evolution models (e.g., 



19971 ). Pure Luminosity Evolution models (in which galaxies form 



at some r edshift and are charac terized by their star-formation history, but no merging oc- 
curs; e.g., 



(.e.g. 



Metcalfe et al. 



Ferguson and McGaugh . 



1996 ), an d models dominated by low-surface-brightness galaxies 

f-ligh t radii that are considerably larger 



1995), all produce ha 



than those observed for 24 27.5. 



Roche et al 



(1998J), who used profile fitting to deter- 



mine half-light radii, also find that at z > 0.35 galaxies are significantly more compact than 
predicted by Pure Luminosity Evolution models (although no evolution in the size was found 
for z^0.35), and more generally, that galaxies at z^2 are more compact than today's L* 
galaxies. 

There is no question that the most attractive explanation for the higher fraction of mor- 
phologically peculiar, and smaller galaxies in the past (beyond z ~ 1 smaller angular sizes 
correspond to smaller physical sizes, essentially irrespective of the cosmological model), is 
in terms of hierarchical galaxy formation and interactions. In the context of this model, the 
highly irregular morphologies are a consequence of interactions (that were more frequent in 
the denser past), direct collisions, and mergers. In this picture, the smaller size objects are 
essentially the "building blocks" of today's galaxies. 

Some attempts have been made to transform this qualitative stat ement into a quantita- 



tive tool. While the Hubble Sequence classification (jHubblel . 11226) has proved extremely 
useful for analyzing the gross morphological properties of nearby galaxies, it becomes rather 
ineffective at high redshifts, when mergers and intense star formation become the rule, rather 
than the exception. 

A method that is quite successful in distinguishing between irregular galaxies forming 
stars stochastically and irregularities induced by mergers is based on color-asymmetry di- 
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agrams (e.g.. 



Bershadv. Jangren. and Conselice 



2000 



Conselicd . 119971 ). In these diagrams, 



galaxies are placed in the \{B — V),A\ plane, where A is an asymmetry parameter, based 
on comparing the galaxy image with its counterpart obtained through a rotation by 180°. 
Formally, the asymmetry parameter is defined by 



.4 



mm 
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E|/ | 



(23) 



where I Q is the intensity distribution in the image pixels, is the intensity distribu- 
tion in the rotated image (by angle 0), and B is the intensity in background pixels 
(with the second term correcting for t he noise). A calibration using nearby galaxies 



(jConselice. Bershadv. and Jangern . 



20q"qT ) shows that all the galaxi.es with A > 0.35 rep- 



resent merging systems. Galaxies with enhanced star formation rates are bluer and may 
become asymmetric due to pockets of star formation. Nevertheless, their asymmetry pa- 
rameter does not exceeed 0.35. 

An examination of the color-asymmetry diagram of the HDF-N shows that the highest 
fraction of mergers is found in the highest redshift range, 1.5 < z < 2.5. Using A > 0.35 
and Mb < — 18 a s criteria for mer gers, gives for the merger fraction as a function of redshift 
/ oc (l + z) 2 - 1±a5 (|Conselicd . 120011 ) . While considerable uncertainties still exist, this suggests 
that the evolution of galaxies may indeed be primarily dominated by mergers. Furthermore, 
the merger fraction appears to be starting to flatten (or possibly even to decline) for z^2. 
As we shall soon see, this may be related to the behavior of the star formation rate with 
redshift. 

Additional information about hierarchical vs. monolithic formation models comes from 
studies of elliptical galaxies. In the former scenario, giant ellipticals form via mergers of 
galaxies of comparable ma ss that had already fprior to merger) us ed up at least some of 
their gas to form stars (e.g.. iKauff mann. White, and Guideroni, 1993 ). In the monolithic see 



iMli 

3g,b 



nario, on the other hand (e.g.. lEggen. Lvnden-Bell. and Sandag' 



0, 



1962; 



Tinslev and Gunn . 



19761 ). ellipticals form at high redshifts via a single collapse (and a concomitant starburst), 



and evolve passively thereafter. The two pure models have rather distinct observational 
predictions. Clearly, in the hierarchical scenario the number density of ellipticals is ex- 
pected to decrease with increasing redshift. This can be contrasted with the predictions of 
the monolithic scenario: the number density of ellipticals should stay fairly constant with 
redshift, but the bolometric luminosity is expected in increase (up to the starburst phase). 
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The obser vations, however, proved to be more ambiguous than one might have hoped. In 



particular, 



Treu and Stiavellil f)1999r ) found that while Pure Luminosity Evolution models in 



which all ellipticals are assumed to form at z ~ 5 over-predict the observed counts, these 
models t end to under-predict the counts when the formation redshift is assumed to be z ~ 2 



fsee also 



Benitez et al. 



1999: 



Franceschini et al 



1998; 



Zepi 119971 ). Thus, the observations 



with both ground-based surveys and HST suggest a picture which is somewhat intermedi- 
ate between the pure hierarchical and monolithic scenarios. Namely, objects that resemble 
elliptical galaxies (red, obeying a i? 1//4 law in their lum i nosity profile) existed already a : 
mode rately high (z£1.5) redshifts (e.g., 



McCarthy et al. 



2001 



Moustakas and Somervilk , 



20021 ). These objects continued, however, to experience mergers until z ~ 1, from whi ch 



Giavalisco 



2002). 



time on the resulting giant ellipticals evolved mostly passively (see also 

The Advanced Camera for Surveys, installed on board HST in March 2002, is expected 
to produce a wealth of high-quality data on galaxy sizes and morphologies. We can therefore 
expect even more significant constraints on hierarchical formation models to emerge in the 
coming months. 

As I noted above, the evolution of galaxies is also intimately related to the way they 
assemble their mass and, concomitantly, to their star formation rate. As it turned out, the 
HDF-N proved to be seminal in the discussion of the cosmic star formation history. 



B. The global cosmic star- formation history 



The rest-frame UV luminosity of galaxies traces nicely their metal production rate, be- 
cause both are produced primarily by massive stars. On the other hand, the conversion 
from a UV luminosity density to an actual cosmic star formation rate is somewhat less 
straightforward, since it involves a knowledge of the Initial Mass Function (IMF). 



Early estimates of the meta^ 
ground-based redshift surveys 
absorption line systems (e.g., 



production rate as a function of red shift, pz{z), relied on 



(Galleeo et al. 



Lanzetta et al. 



1995 



1995 



lillv et al 



Pei and Falll 



I 99I a nd data, from QSO 



19951 ). These estimates 



indicated a monotonic increase in the metal production rate from z = to z ~ 1 (with the 
rate at z ~ 1 being about an order of magnitude higher than at the present). 

One of the seminal results to have come out of the HDF-N was the attempt to esti- 
mate the metal production rate at high redshifts and thereby to generate a continuous plot 
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of the star formation rate as a function of redshift (|Madau et all I1996T ). In the original 
diagram, no corrections were made, for either dust extinction or for surface brightness ef- 
fects. Consequently, the results were presented as lower limits. Both in the original work 
and in some of the subsequent work that followed (some of which included the effects of 
dust extinction), it was found that the star formation rate peaks at a redshift of z ~ 1 



2 anc 



1999: 



decreases or stays nearly con s tant at higher redshifts (e.g ., 



Hopkins. Connolly, and Szalavl . 2000; 



Pei. Fall, and Hauser . 



Calzetti and Heckma 



1999[ 



Steidel et al 



1999 



Fig. 26). The question of whether there truly is a decrease in the star formation rate for 
z^2 has become a focal point of the discussion of the cosmic star formation history. Some 
doubts were raised on the basis of dust extinction and selection effects on one hand, and 
cosmological surface brightness dimming effects on the other. Broadly speaking, the UV 
luminosity density could be underestimated if a significant fraction of the star formation 
occurred in environments obscured by dust or in very low-surface-brightness galaxies. Sev- 
eral at temps were made to correct for dust attenuat ion (using the empirical attenuation 



law of 



Calzetti. Kinney, and Storchi-Bergmann . 



199 



the UV spectral slo pe and the far IR emission (e.g., 



Steidel et al 



), by calibrating the relation between 



Meurer. Heckman. and Calzetti . 



1999; 



1999). Othe rs used constraints obta i ned from extragalactic background radia- 



tion and neutral gas (e.g., 



Calzetti and Heckman 



1999; 



Pei. Fall, and Hauser . 



19991 ). Most 



of these investigations concluded that the star formation rate rises from the present to 
z ~ 1-2 and then stays approximately flat to z ~ 5. 

Another way to address the question of the history of mass assembly in galaxies is to try 
to measure the actual stellar masses of galaxies (as opposed to the rates of star formation; 
ideally one would want to do both). To this goal, observations in the near-infrared are 
typically used, si nce t he near-infrared luminosity traces the stellar mass reasonably well. 

ickinson et al\ ([20031 ) used an infrared-selected sample of galaxies from the HDF-N to 
determine the global stellar mass density, Q*(z), for < z < 3. They found that Q*(z) 
increases with time from z = 3 to the present (Fig. 27). Dickinson et al. concluded that by 
z ~ 1, about 50-75% of the present-day stellar mass density had already formed, but that 
the stellar mass density at z ~ 2.7 was about 17 times lower than today. These observations 
appear to be in clear contradiction with scenarios in which most stars in today's spheroids 
formed at z 3> 2, but the observations are in general agreement with a global star formation 
rate that rises from the present to z ~ 1 and then stays fairly flat. 



41 



A different spin on the cosmic star formation rate has been put bv lLanzetta et al\ ((2002). 
These authors claimed that by neglecting cosmological surface brightness dimming effects, 
previous works have missed a significant fraction of the ultraviolet luminosity density at 
high redshifts. Specifically, since the surface brightness decreases with redshift as (1 + z)~ 3 
(because of the cosmic expansion) , intrinsically faint regions of high- redshift galaxies become 
undetectable. 

Lanzetta et al. designated the unobscured star formation rate intensity (i.e., the intensity 
inferred from the observed rest-frame UV light) by x, and used a distribution function h(x) 
(defined so that h(x)dx is the projected proper area per comoving volume of star formation 
rate intensity in the interval x to x + dx), to estimate the ultraviolet luminosity density 
at high redshifts (including the surface brightness dimming effects). They found that the 
star formation rate density, p s , increases monotonically with redshift to the highest redshifts 
observed (z ~ 8; although in one of the possible corrections for incompleteness p s remains 
fairly flat above z ~ 2). 

The low fraction of stellar mass formed by z = 3 according to the Dickinson et al. 
results appears to contradict evidence for significant star formation occurring at still higher 
redshifts. One possible way out of this conundrum is that the initi al mass function (IMF ) 
at high redshifts is very top-heavy (tilted towards massive stars; iFerguson et all |2002). 
Massive stars dominate the UV luminosity in star forming galaxies, but their contribution 
to the total survi ving (at lower redshifts) stellar mass is relatively low. 

Most recently IStanwav. Bunker, and McMahonl (J2003) used HST's Advanced Camera for 
Surveys, the ground-based Sloan Digital Sky Survey and the Very Large Telescope (VLT) 
to determine the space density of UV-luminous starburst galaxies at z ~ 6. They found a 
lower bound to the integrated, volume- averaged, global star formation rate at z ~ 6, that 
was about six times less than that at z ~ 3-4. The question of the true behavior of p s above 
z ~ 2 remains, therefore, presently somewhat unresolved. 

Fortunately, a more definitive answer may come in the near future, through a combination 
of planned observations with HST and with the Space Infrared Telescope Facility (SIRTF; 
currently scheduled to be launched in April 2003). The Great Observatories Origins Deep 
Survey (GOODS; Principal Investigator M. Dickinson) will produce a very deep image of 
two fields (the HDF-N and the southern deep field observed with the Chandra Observatory) 
with SIRTF at 3.6-24 /im, and will thereby produce a much more complete census of stellar 
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mass at high redshifts. At the same time, observations (of the same fields) with Hubble's 
Advanced Camera for Surveys (Principal Investigator M. Giavalisco; the observations are 
being carried out as these lines are being written) will determine the star formation rates, 
sizes, and morphologies of galaxies. The combined observations will allow for the first time 
for a determination of the evolving mass assembly distribution f(M, M, t) (where M denotes 
the stellar mass and M the star formation rate). Given the fact that under pure luminosity 
evolution, M — J Mdt, a comparison between the observed evolution of the (M, M) phase 
space with time and the evolution obtained from direct integration (of M, to produce M), 
will allow, in principle, for an identification of the role of mergers and interactions. Even 
when the expected observational uncertainties are taken into account, there is no doubt that 
the planned GOODS observations will yield a huge step forward in the understanding of the 
assembly of present-day galaxies, and their morphological evolution (the emergence of the 
"Hubble Sequence"). Furthermore, the planned Hubble Ultra Deep Field with the Advanced 
Camera for Surveys (currently scheduled for July-August 2003) could extend the redshift 

.ose to the tail of the tentative second reionization epoch 
the u niverse (the first reionization having tentatively 



coverage unambiguous! v to z ~ 6, cl< 


(Becker et al. 


2001 


Fan et al 


2001) 



occurred at z 



20+ 10 - 



9 i 



Bennett et al 



2003). Such a study could therefore produce results 
that are not merely incremental in our understanding of the cosmic star formation history, 
galaxy evolution, and the ionization history of the universe, and that can be u sed to place 



meaningful constraints on theoretical models (e.g., 



Somerville and Livio . 



2003). 



VII. THE HUBBLE CONSTANT 



A. A brief background 



Ever since Edwin Hubble's pioneering measurements in the 1920s (and to som e extern 
even before that, since Vesto Slipher's measurement s, that started in 1912; e.g., Slipher, 



1929 



Hubble and Humason . 



1917), we knew that we live in an expanding universe (IHubble . 
19311 1 . In the standard big bang theory the universe expands uniformly, with the recession 
velocity being related to the distance through the Hubble law, v = H^d. More generally, 
based on the "Cosmological Principle" (the assumption that the universe is homogeneous 
and isotropic on large scales), the expansion is governed by the Friedmann equation in the 
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context of general relativity 



r. 



= H 2 



8nGp 



kc 2 Ac 2 

r j + ~ 



(24) 



Here R(t) is the scale factor, H = R/R measures expansion rate (with H , the "Hubble Con- 
stant," giving the rate at present), pu is the mass density, k is the curvature parameter and 
A is Einstein's cosmological constant (which represents the energy density of the vacuum). 
Commonly, the density of matter and that associated with the vacuum are represented by 

the density parameters (at present) Qm = 8ttGpm/3Hq and = Ac 2 /3Hq, with which the 

kc 2 

Friedmann equation can be expressed as — - = Hq(CIm + — !)• The Hubble Constant is 

R o 

thus the key parameter in determining the age of the universe (with VI m and Ha also playing 
a role). Similarly, physical processes such as the growth of structure and the nucleosynthesis 
of light elements (H, D, 3 He, 4 He, Li), as well as critical epochs in the Universe's history, 
such as the transition from a radiation-dominated to a matter-dominated universe, depend 
on the cosmic expansion rate and thereby on the value of Hq. It should therefore come as 
no surprise that the determination of the value of the Hubble Constant became a major 
observational goal for the past eight decades. 



The first value for the Hubble Constant may have actually been derived by 



Lemaitre 



( 19271 ). who, on the basis of Slipher's radial velocity measurements and Hubble's mean 
absolute magni tude for galax ies ("nebulae") obtained Hq = 526 km s _1 Mpc . The next 



set of values bv iHubbld (|l929l ) and Hubble and Humason (who produced a velocity-distance 
relation up to V ~ 20000 km s _1 and obtained Hq = 559 km s _1 Mpc -1 in 1931), were all 
around 500 km s" 1 Mpc -1 , with an uncertain ty stated rath er naively as "of the area of ten 



percent." About twenty years passed before 



Baadd (|1954l ) revised the distance to nearby 



galaxies, recognizing that Hubble confused two classes of "standard candles" (Population I 



Cepheids and Population II W Virginis st ars) , t 
a factor two (a revision suggested also by 



Behr 



lereby reducing the value of H by about 



19511 ). The value of the Hubble Constan t 



first reached the range of values accepted today through the work of Alan iSandaed (J1958). 
Sandage demonstrated that Hubble mistakenly identified H II regions as bright stars, and 
he [Sandage] was able to revise the value further to H ~ 75 km s _1 Mpc -1 (recognizing 
that the uncertainty could still be by a factor 2). 

In the three decades that followed, published values of the Hubble Constant v aried by 
about a factor of two between ~ 100 and 50 km s _1 Mpc -1 . Table 1, adapted from 



Trimble 
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(1997), summarizes the early history of the constant. 

Generally, since redshifts (and therefore radial velocities) can be determined relatively 
readily (this is, of course, not true for the most distant or faintest objects), the problem 
of determining the Hubble Constant has always been a problem of determining accurate 
astronomical distances. The availability of new instrumentation, and the Hubble Space 
Telescope in particular, have allowed for a dramatic improvement in distance determinations. 

B. Distance indicators and methods 

Direct trigonometric parallaxes that use the Earth's orbit around the Sun as a baseline 
for triangulation can only be used to the nearest stars. Consequently, othe r distance in 
dicatqrs and method s had to be used for extragalactic distances (see, e.g. 



1992; 



Trimble 



Jacobv et al. 



1997L for excellent reviews). The most common of these employ "standard 
candles," geometrical properties, physical properties, or various correlations to determine 
distances. 

Standard candles are simply based on the fact that the flux of radiation decreases as 
an inverse square law. Objects with either a constant luminosity or whose luminosity can 
be related to a distance-independent measurable property (such as an oscillation period) 
are good standard- candle candidates. The best known and probably most reliable in this 
class are the Cepheid variables. Their potential as standard candles on the basis of their 
period-lumin osity (P-L) rela tion was first recognized by Henrietta Leavitt in 1912, and they 



were used by iHubbld (|1925[ ) to determine distances to Local Group galaxies. The physical 
processes responsible for the P-L relation are broadly understood. Near ionization zones 
(in this case, primarily He + 1± He ++ ), gas can absorb heat under compression and release it 
after maximum compression. This leads to an instability strip in the effective temperature- 
luminosity plane, which, in the case of Cepheids, is very narrow in temperature. The 
pulsation period depends on the mass and radius as 

1 



~ M~2R2 . (25) 

The luminosity of the star (which is determined ultimately by nuclear reactions that depend 
on the density and temperature) is proportioned to a power of the mass, L ~ M k . Since, 
however, we also have (for black body radiation) L ~ R 2 T 4 , we obtain P ~ L~*~T~ 3 , or a 
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period-luminosity-color relation. One of the "Key Projects" of HST has been to measure Ho 
based on a Cepheid calibration of a number of secondary distance determination methods. 
The results of this project have been described in a series of some 30 papers (see 



Freedman et al. 



2001 



and references therein). The main goals of the project have been: 
(i) To discover Cepheids in a sample of relatively nearby galaxies (with distances £j 20 Mpc) 
and to determine distances to these galaxies, (ii) To determine Ho through several secondary 
distance indicators, to all of which Cepheid calibration is applied; and perhaps most impor- 
tantly, (iii) To determine the uncertainties in all the methods by comparing the distances 
obtained from them to Cepheid-based distances. As a part of the project, the uncertainties 
in the Cepheid P-L relation itself (and its dependence on other factors such as metallicity) 
have been investigated. 

The Key Project used Cepheid calibration to 31 galaxies, of which 18 have been observed 
and analyzed in the context of the project, and to which archival data, and data on the 
nearby galaxies M31, M33, IC 1613, NGC 300, and NGC 2403 have been added. 

The Key Project used the following secondary methods based on Cepheid distances: 
Type la Supernovae, the Tully-Fisher Relation, the Fundamental Plane for elliptical galaxies, 
Surface Brightness Fluctuations and Type II Supernovae. Let me describe very briefly the 
physical basis for each one of these methods. 



1. Type la supernovae 



Type la supernovae at peak brightness are extremely bright, with Mb 



Mi 



v 



— 19. 3 + 5 log(ifo/60), and they show a relatively low dispersion, <t(Mb) ~ 0.33 fe.g.. lBranch . 
19981 ). Furthermore, there exists a relatively tight correlation between their peak lumi- 
nosity and light- cury e shape for rat e of d e cline, with b righte r supernovae declining more 



slowly; 



Hamuv et al 



1996 



Phillip, H 



Riess et al, 



1996). The homogeneity may be 



related to the fact that Type la supernovae represent thermonuclear disrupti ons o f mass - 



accreting white dwarfs, when the latter reach the Chandrasekhar limit (e.g. 



Livio, 



Arnett , 



2001). 



2001) 



The luminosity-light-curve relation may be the result of the following (e.g., 
The peak luminosity of a supernova Type la is proportional to the mass of 56 Ni that is 
produced. A higher mass of 56 Ni, however, also results in more heating and concomitantly 
a higher opacity (due mainly to UV lines). Consequently, a slower development of the light 
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curve results. 



2. The Tully-Fisher relation 



For spiral galaxies, th at are kuown to have flat rotation curves, there is an observationally- 
determined relationship ([Tullv and Fisherl . Il977l ) between the total luminosity and the max- 
imum rotational velocity (both corrected for inclination effects), of the fo rm L ~ V 3 ^ (in 



the I band). The scatter around this relation is about ±0.3 mag (e.g., 
19971) . 



Giovanelli e.t al. 



ation. Very broadly , 



Eisenstein and Loeb . 



There is no precise physical understanding of the Tully-Fisher re 
a sim ilar relation can be obtained from the following argument (e.g., 
19961 ). Consider a galaxy (of mass M) collapsing from a spherical cloud. At epoch £ co u, the 
turnaround radius R t is 



Rt 



After virialization, the energy is 



E ~ Ma 2 ~ GM 2 /R t 



(26) 



(27) 



where a is the velocity dispersion. Combining the above and assuming that the galaxy forms 
from a single collapse gives 

a~(M/t coll )i . (28) 

Therefore, if all galaxies collapse at the same epoch and the mass-to-light ratios, M/L, do 
not vary significantly, we obtain L ~ a 3 . 



3. The fundamental plane for elliptical galaxies 



Large spectrophotometri 


c survevs 


Dioreovski and Davis 




1987 




Dressier et al. 



(e.g., 



19871 ) revealed that, for elliptical galaxies, 



a tight correlation exists between the effective radius, R e , the effective surface brightness, 
SB e , and the central velocity dispersion a, of the form (the "fundamental plane") 



logi^ = a loga + (3SB e + 5 



(29) 
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Here R e is in kpc, a in km s -1 , and SB e in mag arcsec" 2 . The value of 5 depends on Ho, 
since the calculation of the effective radius in kpc uses the Hubble constant. 

The physical origin of the fundamental plane relation can be understood on the basis 



Treu et al 



of the following simple considerations (e.g., 
(virial) mass by 

M = a 2 R e /G . 

Let us also assume that the mass-to-light ratio satisfies 

M 



We then obtain 



M 6 



20011 ). We can define an effective 

(30) 

(31) 



(32) 



which reduces to the fundamental plane relation for reasonable values of S (8 ~ 0.25). 



4- Surface brightness fluctuations 



The me thod of surface brightnes s fluctuations was developed by 



1988 


) and 


Tonrv et al. 


(1997. 


2000) 



Tonrv and Schneider 



2000). The method basically makes use of the obvious fact 
that the ability to resolve stars within galaxies is distance dependent. More specifically, for 
every region of a galaxy one can measure the average flux per pixel, g, and the pixel-to-pixel 
rms, a. Since the flux obtained in a pixel is received from N stars of average flux /, we have 
g = Nf and a = f. A galaxy which is twice as distant appears twice as smooth as the 
closer galaxy. Consequently, the average stellar flux is given by / = cr 2 /g (and / scales as 
the inverse of the square of the distance). 



5. The Expanding Photosphere Method of Type II supernovae 

Type II supernovae result from the collapse of stars more massive than about 8 M Q . 
Generally, Type II supernovae are fainter than Type la supernovae, and they also exhibit 
a considerably larger range in their luminosities, making them poorer standard candles. 
Nevertheless, Type II supernovae have been used as distance indicators, through an applica- 



tion o f the Expanding Photosphere Method fe.g.. lKirshner and Kwan . 
19941) . 



1974 : 



Schmidt et al. 
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The basic idea is simple. The angular size of the photosphere is given by (for z<l) 

where T is the color temperature, f\ is the flux density, B\(T) is the Planck function, and (\ 
represents the dilution effects of scattering atmospheres (derived from model atmospheres). 
The photospheric radius is given by 

R = v{t-t ) + R , (34) 

where v is the expansion velocity (measured from the absorption minima of optically thin 
lines), and the initial radius, R Q , can be neglected at all but the earliest times. Combining 
the above yields 

t = D(-)+t , (35) 



v . 

making it possible to determine both the distance and the time of the explosion from a few 
measurements of t, v, O. 



C. The results 



Table 2, adapted from 



Freedman et al. 



( 20011 ). lists the values of H obtained from the 



different methods, based on the Cepheid distances of the Key Project. Of the five calibrated 
methods, it is clear that the Fundamental Plane is somewhat of an outlier. Combining these 
results, the H Key Project obtained the value H Q = 72 ± 3 ± 7 km s _1 Mpc -1 , where the 
first quoted error is random and the second is systematic. Using three different weighting 
schemes, all the results were found to be consistent with Ho = 72 ± 8 km s _1 Mpc -1 . 

In an independent work, Allan Sandage, Gustav Tammann, Abijhit Saha, and collabora- 
tors used a Cepheid calibration of the peak brightness of Type la Supernovae (using HST), 
and thereby determined Hq directly from the Hubble diagram of the latter. This effort 
resulted to date in nine supernovae la with "normal" spectra (i.e., the spectra characteriz- 
ing about 60 % of all Type la s upernovae) to which Cepheid distances are known (Table 3, 



adapted from 



Saha et al. 



20011 1. After corrections based on the decline rate and colors, the 



weighted average of the absolute magnitudes of the calibrators was fitted to a sample of 35 
more distant, "normal" Type la supernovae. After further correcting for some systematic 
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errors. 



Saha et al 



( 20011 ) obtained Hq = 58.7±6.3(internal) km s _1 Mpc -1 (for cosmological 
parameters Qm = 0.3, Q\ = 0.7). 

Unlike the uncertainty by a factor of two that has plagued this field for decades, therefore, 
the observations with HST have reduced the uncertainty in the value of the Hubble constant 
to about 15%. 

The values obtained by both the Key Project and the Sandage-Tammann-Saha team are 
also consistent with other recent measurements based on combining the Sunyaev-Zeldovich 
effect with x-ray flux measurements of clusters, and on time delays in gravitational lensing. 
The systematics in these methods are still estimated to be at the 20- 25% level, and the 
value obtained for th e Hubble constant is H ~ 60 km s _1 Mpc -1 (e.g., 



2000; 



Schecter, 



200(4) 



Reese et al 



2002, 



Analysis of Sunyaev-Zeldovich effect ([Sunvaev and Zeldovich 



1970) and x-ray data pro- 



vides a direct method for determining distances to galaxy clusters. Clusters of galaxies are 
known to contain hot intracluster gas (at kT ~ 10 keV), trapped in the clusters' potential 
wells. Photons from the cosmic microwave background that pass through a cluster have a 
finite probability (optical depth r ~ 0.01) to interact with energetic electrons in the intra- 
cluster gas. The inverse Compton scattering that ensues boosts the energy of the microwave 
background photon, generating a small distortion (a decrement in frequencies ,$,218 GHz 
and an increment above this value) in the spectrum of the microwave background. The 
Sunyaev-Zeldovich effect is proportional to the integral of the pressure along the line of 
sight (/ n e T e dl). Since the x-ray emission from the intracluster medium is proportional to 
a different power of the density, / n^Adl (where A is the cooling function), a combination 
of the two measurements (given some assumptions about the cluster geometry) can be used 
to determine the distance to the cluster, independently of the distance scale ladder that is 
based on standard candles. 

Finally, the Wilkinson Microwave Anisotropy Probe (WMAP) found a Hubble constant 
of Hq = 72 ± 5 km s -1 Mpc" 1 . When the WMAP data were combined with the Key 
Project results, finer-scale cosmic microwave background experiments, and other large-scale 
structure and Lyma n a forest data, the best-fit value of the Hubble constant was H = 
711% km s" 1 Mpc" 1 jspererel et allhooA . 

Type la supernovae have played a key role not only in the determination of the age of 
the universe (through Hq), but also in the determination of the universe's geometry, and the 
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dynamics of the cosmic expansion. This came about through a combination of ground-based 
and HST observations in which both careful planning and serendipity played a part. 



VIII. THE ACCELERATING UNIVERSE 



In 1998, two teams of astronomers, wo rking independen t 



expansion of the universe is accelerating flPerlmutter et al 



y, presented evidence that the 



1999; 



Riess et al 



19981) . This 



evidence was based primarily on the faintness (by about ~ 0.25 mag) of distant Type la 
supernovae (SNe la), compared to their expected brightness in a universe decelerating under 
its own gravity. 

The first suggestions that high-redshift Type la supernovae could be used to determine the 



rate o f cosmic deceleration came in the late 1970s (jColgate . 



1979; 



Tammann 



1979; 



Wagoner, 



19771 ) , but the actua! 



discovery of an even moderate redshift (z = 0.31) SN la (SN 1988U) had 



to wait for a decade (jNorgaard-Nielsen et a l. 1989). The samples of high-redshift supernovae 



started to become sufficiently large to place constraints on cosmological parameters through 
the efforts of the Supernova Cosmology Project, led by Saul Perlmutter, and the High-z 
Supernova Search Team, led by Brian Schmidt. By 1998 the two teams gathered sufficient 
data to be able to show that the universe is characterized by a m atter density paremater 
satisfying fl M < 1 (i.e. , that the universe is not closed by matter; 



Perlmutter et al 



Garnavich et al 



1998; 



19981 ). However, the real shocker was still to come. 



A. The observations 



The two supernova search teams use a similar method to find their supernovae. They 
take two deep images separated by about a month, subtract the first-epoch image from the 
second-epoch one, and search for sources above a certain threshold in the difference image. 
Once a candidate SN is identified, the SN type is determined by its spectrum (if that can be 
taken; in a few cases one has to rely on the host galaxy type — only SNe la were found so far 
in ellipticals). Supernovae at relatively high redshift are then monitored photometrically to 
construct their light curves. An image of the host galaxy is obtained at a later time (after 
a year or more), and subtracted to obtain an accurate measurement of the SN brightness. 

The Hubble Space Telescope has proved to be crucial especially for the highest-redshift 
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supernovae. There, the ability to resolve and pinpoint the SN location on the host (includ- 
ing in cases in which the SN was found close to the galactic nucleus) was essential for a 
correct determination of the SN magnitudes. With samples of a few dozen SNe la in hand, 
the two teams compared their measured distances (derived from the luminosity-distance 
relation, F = L/AttDI; where Dl is the distance and L, F are the intrinsic luminosity 



and observed flux, respectively) wi th the distances expected 
different cosmological models fe.g.. lCarroll. Press, and Turner! . 



or the observed redshifts, for 



1992). The latter is given by 



D L = cHo l (l + z)\l-Q M - fi A r 1/2 sinn {|1 - Q M - Ov| 1/2 x 

J Z dz[(l + zf{\ + Q M z) - z(2 + z)Q A ]- 1/2 } (36) 

where sinn denotes sinh for Qm + < 1 and sin for Qm + > 1- The results for the 
likelihood of the cosmological parameters Qm and Q\ are shown in Fig. (28). As can be seen 
from the figure, the results favor values of Om — 0.3, Qa — 0.7 and a negative deceleration 
parameter (a is the scale factor, fl is the sum of today's energy densities, and w = Pa/ Pa 
characterizes the dark energy "equation of state"; see section VIII D) 

_ -a(to)a(t ) Q 3 
q °= a*(t ) = T + 2^ A< ° ' (3?) 

corresponding to an accelerating universe. 



B. Alternative interpretations 



Clear ly the interpretation of an accelerating cosmic expansion cannot go unchallenged 



(see also 



Riess 



2000; 



Turner . 



2000). The two most serious challenges to acceleration are, in 



principle, evolution/progenitors and dust extinction. 

The possibility that SNe la undergo some type of cosmic evolution cannot be easily dis- 
missed. After all, almost everything else (e.g., the host galaxies themselves, the metallicity) 
evolves. The inference of acceleration, on the other hand, assumes that the local and high- 
redshift supernovae are drawn from the same statistical sample. I should explain that by 
evolution, I do not mean that somehow the processes governing supernova explosions are 
necessarily changing with redshift, but rather that when observing the high-redshift uni- 
verse we are sampling younger galaxies, and consequently the population from which the 
supernovae are drawn is typically younger. 
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Since SNe la represen t thermonuclear disruptions of mass-accreting white dwarfs 



f Hil lebrandt and Niemever , 



2000; 



Livio 



2001), changes in the C/O ratio in the white dwarfs 



(as a result of lower metallicities in t he distant pas t ), for example, could produce some in 



homogeneity in SNe la light curves ((Umeda et al 



19991 ) . However, observations tend to 



show that even if such an effect exists, it is insignificant. In fact, existing o 



local galaxies (e.g.. ICappellaro et al. 



Dservations of 



1997 



Riess et al. 



1999; 



van den Ber gh , 



19941) already 



span a wider range of metallicities and host properties than that expected on the average 
between the sample of galaxies at z ~ and z ~ 0.5. 

To my knowledge, there is only one known evolutionary effect that is physically meaning- 
ful, that can mimic accelerated expansion. This is the effect of the metallicity on the density 
at the point of central carbon ignition (the trigger of SNe la). Generally, a lower m etallicity 
(as expected at high-z) will result in a lower central density (|Nomoto et a/I Il997l ). This is 
because a lower metallicity results in a lower abundance of 21 Ne, that is responsible for much 
of the neutrino cooling (via the so-called local URCA shell process, involving the 21 Ne- 21 F 
pair). A lower metallicity therefore reduces the cooling and leads to an earlier ignition. Due 
to the lower white dwarf binding energy, the light curve exhibits a more rapid development, 
and a lower implied maximum brightness. However, as I noted above, the local (low-z) 
sample of galaxies does not appear to exhibit any significant metallicity-dependent effect. 

Dust extinction is another natural candidate for making the distant supernovae appear 
dimmer. In fact, the observed decrease in the brightness (~ 0.25 mag) of the distant sample 
would require only a ~ 25% increase in the extinction. Both teams used colors to correct 
for dust extinction (like in sunsets, ordinary dust also reddens the light). Furthermore, 
observations from the optical to the infrared of SN 1 999Q (at z = 0.4 6) showed that a 



20001 ) , making the dust 



large extinction to this supernova is highly improbable ((Riess et al 
hypothesis untenable. 

Another potential uncertainty is associated wi th the progen itors of SNe la. The progenitor 
systems of SNe la are not known without doubt (|Liviol . 120011 ). In particular, the two leading 
progenitor -system candidates are double whi t e dwa rf systems (that merge to produce the 



explosion; 



Iben and Tutukov 



1984; 



accretes from a normal companion ( Nomot 



Webbink, 



1984) and systems in which the white dwarf 



Whelan and Iben , 1973 ) . The possibility 



therefore exists, in principle, that the local and distant samples are dominated by different 
progenitors, with one class producing somewhat dimmer supernovae. This possibility should 
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certainly be considered, especially in view of the fact that 



Li et al. 



( 20011 ) find a relatively 



hi gh (~ 40%) fraction o f pecu liar SNe la among the local sample. 

Yungelson and Livid ([20001 ) have shown that while it is possible, in principle, that one 
class of progenitors (e.g., double white dwarf systems) dominates the local sample and 
another (a white dwarf with a normal companion) the high- z one, such a transition is not 
likely, because of the following reason. If such a transition were to occur, one would expect 
that around the transition point (at z ~ 1), SNe la would be composed of an equal mix of 
the two classes. This is inconsistent with the observ ations, that show that the high-z sample 
is actually extremely homogeneous (jLi et all 120011 ). Consequently, it is highly unlikely that 
the interpretati on of an acce lerating universe is an artifact of the existence of two classes of 



progenitors (see 



Livio 



2001 



for a more complete discussion). 
Other effects, such as those resulting from gravitational lensing, will require some further 
consideration once larger databases of high-redshift supernovae become available. Generally, 
it can be expected that most lines of sight to detected supernovae do not pass through large 
mass concentrations. Consequently, light paths are being bent out of the line of sight- 
resulting in deamplification. In rare cases however, the line of sight may cross significant 
matter concentrations resulting in strong amplification. The picture that emerges, therefore, 
is that of a brightness distributio n in which the peak is shifted toward a lower value, but 



with a tail of bright objects (e.g., 



Holz . 



1998). Effects of this type can be averaged out once 



deep data from many lines of sight become available, perhaps with the proposed wide-field 
imager, the Super Nova/Acceleration Probe (SNAP). 



C. Supernova 1997ff serendipity and careful planning 



In the spirit of "extraordinary claims require extraordinary proof," the magnitude of the 
discovery of an accelerating universe requires an extremely careful elimination of potential 
contaminating astrophysical effects. A pervasive screen of "gray" dust, for example, that can 
dim the light while leaving little imprint on the spectral energy distribution, has been sug- 



geste d as an alternative explanation for the unexpected faintness of high-;? SNe la ([Aguirre , 



19991 ) . While gray dust, e.g., quasispherical grains larger then 0.1 /xm, could be made to 
reproduce some of the observations, some measuremen ts do disfavor a 30% visual opacity 



due to gray dust at the ~ 2.5a level (jRiess et al. 



2000). Nevertheless, both the possibility 
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of gray dust and, in particular, the ever-existing challenge of luminosity evolution, did leave 
the case for an accelerating universe somewhat short of compelling by 2001. The strongest 
evidence so far against alternatives to an accelerating universe came from SN 1997ff. 



In 1997, Ron Gilliland and Mark Phillips reobser yed the Hubble Deep Field Nor t 



i with 



HST, with the goal of detecting high- z supernovae f Gilliland. Nugent and Phillip; 



1999). 



They discovered two supernovae: SN 1997fg, at a redshift of z = 0.95 was found in a 
late-type galaxy, while SN 1997ff was found to be hosted by an elliptical galaxy. The ini- 
tial, photometrically determined redshift, placed the latter galaxy at z ~ 0.95-1.32 with 
co nsiderable uncert a inty. Serendipitously the Guaranteed Time Observer (GTO) program 



of 



Thompson et al\ ((1999) included imaging the host of SN 1977ff with the near infrared 



camera, NICMOS, on board HST. One of these exposures, amazingly enough, was actually 



taken wit 
program 



lin hours of the disco very of the SN. Furthermore, six months after the GTO 



Dickinson et al 



(2003J, in preparation; GO 7817) reobserved the field with NIC- 
MOS. The combination of all the HST and ground-based observations (in the U, B, V, I, 
J, H, and K bands) all owed for an improved photometric redshift determination for the 



host of z = 1.65 ± 0.15 ( Budavari et al 



20001 ). making SN 1997ff the farthest known su- 



pernova to date. A probability density function based on the partial supernovae light curve 
in three bands also gave z ~ 1.7, as did a tentative determina tion of the redshift based on 



Riess et al. 



2001 



for a full 



the spectrum of the host taken with the Keck telescopes (see 
description). 

The importance of detecting a SN la at such a high redshift cannot be overemphasized. 
Even with the suggested value of the vacuum energy density of Q\ ~ 0.7, at z ^ 1 attractive 
gravity would still have dominated over the cosmic repulsion, resulting in a decelerated ex- 
pansion. Consequently, SNe la in the redshift range z ~ 1-2 should appear brighter relative 
to SNe in a coasting universe. On the other hand, any effect that increases monotonically 
with redshift, as would be expected from dust extinction and simple evolutionary effects, 
would predict that the SNe la would be dimmer with increasing redshift. SN 1997ff was 
found to be brighter by ~ 1.1 mag (at the > 99.99% confidence level) than expected for an 
alternative source of dimming (e.g., dust or evolution) beyond z ~ 0.5. Figure 29 shows a 



redsh ift- dist ance relation in w hich the points are redshift-binned data from 



(119991) and 



Riess et al. 



Perlmutter et al. 



(1998), together with a family of curves for flat cosmological models. 



The transition from a decelerating to an accelerating phase occurs (if the "dark energy" is 
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represented by a cosmological constant; see Section VIII D) at z tr = — 1. As 

can be seen from the figure, the only cosmological model that is consistent with all the data 
(including SN 1997ff) is one in which tt M ~ 0.35, fi A ~ 0.65. 

Measurements of the power spec t rum of the co s mic microwave background (e.g., 



Abroe et al. 



2002; 



de Bernardis et al 



2002; 



Hu et al. 



2001 



Netterfield et al. 



2002) pro- 



vided strong evidence for a geometrically flat (Om + — 1) universe. When combined 
with several estimates of Qm from mass to light ratios, x-ray temperature of i ntracluster gas 
numbers, and dynamics of clusters of galaxies, all giving Om — 0.2-0.3 (e.g ., Bahc all et al. 



200 



0; 



Ba hcall. Fan, and Cen 



1997 



Carlberg et al 



1996; 



Strauss and Willick . 1995), the in- 



escapable conclusion was that there is a dark energy component of Q\ ~ 0.7 — consistent with 
the value found from high-redshift supernovae. Most recently, the best fit cosmological model 
has been determined by the Wilkinson Microwave Anisotropy Probe (WMAP). In combi- 
nation wit h other large-scale structure work, the W MAP results imply f2 to t = 1-02 ± 0.02, 



Cl A ~ 0.73 (^Bennett et al 



2003; 



Spergel et al 



2003). 



D. The nature of "dark energy" 



It is beyond the scope of the present a rticle to discuss dark ene rgy in detail and the 



reader is referred to the excellent review by 



Peebles and Ratra 



(2002). However, since this 



is arguably the most dramatic discovery involving HST I would like to make a few points. 
In general relativity, the stress-energy tensor of the vacuum, T^ v , can be written as 



T, 



(38) 



where p\ is constant (proportional to Einstein's cosmological constant) and g^ v is the metric. 
In a Minkowski flat spacetime this can be written simply as the "equation of state" 



Pa = -Pa • 



(39) 



More generally, if the equation of state is wr i tten as (for example, j Canuto et all Il977 



Ratra and Peeble; 



1988; 



Sahni and Starobinskv 



2000; 



Turner and Whit 



19971 ) 



P A = wp A 



(40) 



then the dark energy density behaves (for a constant w) like p\ ~ a 3 ( 1+u ') (where a is the 
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scale factor). Since the scale factor satisfies the equation (in units in which c = 1) 

! = ^(P + 3P) , (41) 
a 6 

the dark energy contribution will result in an accelerating universe if w < — |. 

It is interesting to note that from a purely fluid dynamical point of view, any equation 
of state of the type P = wp is in fact dynamically unstable for any negative value of w 
except for w — — 1 (corresponding to a cosmological constant). This can be easily seen from 
perturbations on the equations for momentum and energy conservation (< p >, < P > 
denote mean values) 

5p+ (< p > + < P >)V-v = (42) 
(< p> + < P >)$+ \^P\V5p = . (43) 



dP 

Defining the speed of sound, c\ = — — , these two equations can be combined to form 

dp 

5p = c 2 s V 2 5p , (44) 

which is clearly unstable for negative c 2 s (except for w = —1, when the combination < p > 
+ < P > vanishes, as expected from the fact that the vacuum is the same for all inertial 
observers) . 

Thus, formally speaking, from a fluid-dynamical point of view only a cosmological con- 
stant yields a stable solution. 

One can certainly take the dark energy to be associated with a uniform scalar field <fi (a 
"quintessence" field; see below), that rolls down a potential V((f>) at a rate determined by 



+ 3#0 + V\<t>) = , (45) 



dV 

where V'(0) = —r- . In this case, 
d<p 



P<t> = \tf + V{4>) (46) 

For a sufficiently slowly varying dynamical component, such that the kinetic energy is much 
smaller than the potential, 2 V(0), one even obtains a field energy that mimics the 
cosmological constant (P^ ~ —p<j>)- Generally, quintessence solutions allow for different 
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equations of state, and even values of w that are time variable. Nevertheless, I find the sta- 
bility property expressed in eq. (44) sufficiently intriguing to warrant a deeper examination 
of the cosmological constant possibility, before that is abandoned in favor of other dynamical 
solutions. This point of vie w received some furth er support from the recent determination 



by WMAP that w < -0.78 ((Bennett et al 



200.4 



There are two main problems with the implied energy density pv of the dark energy (e.g., 



Weinberg 



20011 ). (i) Why is its value not ~ 120 orders of magnitude larger (as expected 



from fluctuations in the gravitational field up to the Planck scale)? (ii) Why now? (Namely, 
why Q\ ~ Qm now, even though Q\ may be associated with a cosmological constant, while 
VLm declined continuously from the initial singularity to its present value). 

An interesting curiosity to note is that even though taking graviton energies up to Planck 
scale, Mp, misses the value of the dark energy density by ~ 10 120 , and taking them up to 
the supersymmetry scale, M SU sy, misses by ~ 10 55 , a scale of My ~ (M SU sy/^p)^susy 
actually does give the right order of magnitude! While I am not a ware at present of any 



Arkani-Hamed et al. 



theor y that produces this scale in a natural way (although see, e.g 
20001 ). I find this coincidence worth following up. 

Much of the efforts to resolve the above two problems revolved not around a cosmological 
constant, but rather around the behavior of quintessence fields. In particular, attention 
has concentrated on "tracker" solutions, in which the final value of t he quintessence energy 
density is independent of fine-tuning of t he initial conditions (e.g 



2000; 



Albrecht and Skordis 



Zlatev. Wang, and Steinhardtl . ll998r ). For example if one takes a potential of the form 



V(cj)) = <p- a M 



(47) 



where a > and M is an adjustable constant, and the field is initially much smaller than 
the Planck mass, then for pu 3> V(</>), 2 , we find that pM decreases initially faster (~ t~ 2 ) 
than 0(t)(~ t-2«/(2+«)). 

Eventually, however, a transition to a /^-dominated universe occurs (and p^ decreases as 
£- 2 /( 4 + Q )) # i n other words, the quintessence answer to the question: why is the dark energy 
density so small? is simply: because the universe is very old. Nevertheless, simple potentials 
of the form (47) do not offer a clear solution to the "why now?" problem. In fact, to have 
P<t> ~ Pm (and of order of the critical density) at the present time requires fine-tuning the 
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parameters so that (for example. IWeinberd . 12001) 

M 4+Q ~ (ZTiG)- 1 -* 12 ^ 



(48) 



with no simple explanation as to why this equality should hold. 

In order to overcome this requirement for fine-tuning, some versions of the quintessence 
models choose potentials in which the universe has periodically been accelerating in the past 
(the dark energy has periodically domin ated the energy density; e.g., Dodelson, Koplinghat 



& Stewart 2000). 
the form 



Dodelson et al 



(2000J) have shown, for example, that with a potential of 



V(4>) = V a exp (-\(j)V8^G) [l + A sin (z/0^8ttG 



(49) 



one can obtain solutions in which the dark energy density tracks the ambient energy density 
of the universe and satisfies observational constraints. 



A different approach to the proble ms associated with the dar 



throu gh anthropic considerations (e.g., 



Kallosh and Linde 



2002; 



t ene rgy 



Vilenkin 



density has been 



1995 



Weinberg . 



20011 ). The key assumption in this class of models is that some constants of nature, in par- 
ticular the cosmological constant, and possibly the density contrast at the time of recombi- 
nation, cr rcc , are in fact random variables, whose range of values and a priori probabilities are 
nevertheless determined by the laws of physics. In this picture, some values of the constants 
which are allowed in principle, may be incompatible with the very existence of observers. 
Assuming a principle of mediocrity — that we should expect t o find ourselves in a universe 



typica l of those that allow the emergence of intelligent life — iGarriga. Livio. and Vilenkin 



(2000J) were able to show that the "why so small?" and "why now?" questions find a nat- 
ural explanation. While I personally believe that anthropic considerations should only be 
used as a last resort, there is no denial of the fact that most versions of "eternal inflation," 



Guth. 


1981; 


Linde . 


1986; 


Steinhardt . 


1983; 


Vilenkin . 


1983) 



an infinite ensemble is believed to exist, th e problem of defining probabilities on it, and the 
conce pt of "mediocrity" become real (e.g., 
1998)). 



Linde. Linde. and Mezhlumian 



1995 



Vilenkin . 



In spite of the apparent "success" of the anthropic argumentation in solving the two main 
problems associated with dark energy, the search for a fundamental explanation is not, and 
should not, be abandoned. An interesting, if speculative, new direction is provided by 
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alternative theories of gravity. For example, models in which ordinary particles are localized 
on a three-dimensional surface (3-brane) embe dded in infinite volume extradimension s to 
which gravity can spread, have been developed ()Deffavet. Dvali. and Gabadadzd . 12002). In 
a particular version of these models, the Friedmann equation (eq. 24) is replaced by 

2 



H 2 



k 



P 



1 

4r? 



(50) 



\ 3M| 4r c 2 2r c 2 

where e = ±1, and r c represents a crossover scale. At distances shorter than r c (which can be 
of astronomical size, e.g., r c ~ cHq 1 ), observers on the brane discover the familiar Newtonian 
gravity. At large cosmological distances, however, the force-law of gravity becomes five- 
dimensional (as gravity spreads into the extra dimensions) and weaker. The dynamics are 
governed by whether p/M 2 is larger or smaller than 1/r 2 . While highly speculative at this 
stage, some aspects of these alt ernative theories of gravity can be experim entally tested (e.g., 



by lunar ranging experiments; 



Dvali. Gruzinov. and Zaldarriaga . 



2002). Also, the relative 



lack of power on large scales (in particular, absence of any correlated s ignal on scales large r 
than 60 degrees) found by the Wilkinson Microwave Anisotropy Probe ([Spergel et all 12003). 
may indicate a breakdown of canonical gravity on large scales (or, more speculatively, a finite 
universe) . 

Luckily, observations that are already being performed with HST, and observations with 
the proposed Super Nova/ Acceleration Probe (SNAP) will provide some crucial information 
on the nature of dark energy. In particular, if the equation of state parameter w is constant, 
then a transition from a decelerating phase to an accelerating phase occurs at 



[(l + 3 W )(fi M -l)/fi M ]- 1/3w -l 



(51) 



For example, for Qm = 0.27 (jBennett et all I2003T ) and a cosmological constant (w = — 1), 



this gives z tr ~ 0.76. Equation (51) shows that if f2 M can be determined independently 
to within a few percent (a goal that can perhaps be achieved by a conbination of cosmic 
microwave background and large scale structure measurements), then observations of a large 
sample of Type la supernovae in the redshift range 0.8-2 can place very meaningful con- 
straints on the dark energy equation of state. The situation is more ambiguous if w is 
time-dependent, but even then, a sufficiently large sample of SNe la will give useful infor- 
mation. At the time of this writing, a search for SNe la conducted in the GOODS fields 
has already yielded 10 SNe la, with six of them in the redshift range z ~ 1-2 (Principal 
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Investigator: Adam Riess). This, together with a few current ground-based searches at lower 
redshifts, is an excellent first step. The survey proposed by SNAP, of 15 square degrees, 
could yield light curves and spectra for over 2000 SNe la at redshifts up to z ~ 1.7. With 
this sample in hand, one could reach the elimination or at least control of the sources of 
systematic uncertainties needed to discriminate between different dark-energy models. 

IX. EPILOGUE 

The Hubble Space Telescope, together with many other space-based and ground-based 
observatories, have given us a much clearer picture than we ever had before of the universe 
we live in. We now know not only that other planetary systems are quite common, but that 
the cosmic expansion is accelerating. We know the age of the universe, its geometry, and its 
composition. Observations of stars start to reach the level of detail previously characterizing 
only solar physics. A general picture of the formation of structure in the universe, and many 
clues for how galaxies form and evolve, start to emerge. We know that supermassive black 
holes reside at the centers of most galaxies. This is, however, far from marking the end 
of HST's contributions. The recent installation of the Advanced Camera for Surveys, and 
the two future (currently planned for 2004) instruments, the Cosmic Origins Spectrograph 
(COS) and the Wide Field Camera 3 (equipped with an infrared channel), promise that the 
coming seven years will be at least as exciting as the first thirteen were. 
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TABLE I The early history of the Hubble constant. Unless otherwise shown, stated error bars 
were 10% or less, or not given. 



When 


Who 


Numerical value 


1927 


Lemaitre 


600 


1929 


Hubble 


530-513-465 


1931 


Hubble & Humason 


558 


1936 


Hubble 


526 


1946 


Mineur 


330 


1951 


Behr 


250 


1952 


Baade, Thackeray 


270 


1956 


Humas, Mayall & Sanda^ 


re 180 


1958 


Holmberg 


134 


1958 


Sandage 


150-75-38 


1959 


McVittie 


227-143 


1960 


Sersic 


125 


1960 


van den Bergh 


125 


1960 


van den Bergh 


125 


1961 


Ambartsumyan 


140-60 


1961 


bandage 


113-85 


1964 


de Vaucouleurs 


125 


1968-69 


de Vaucouleurs 


100 


1969 


van den Bergh 


110-83 


1968-76 


Sandage & Tammann 


50 


1972 


Sandage 


55 


1979 


de Vaucoulers 


100 
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TABLE II Hq from secondary methods (the Key Project) 



Error 
(random, systematic) 
Method H (%) 

36 Type la SN, 4000 < cz < 30, 000 km s" 1 71 ±2 ± 6 

21 TF clusters, 1000 < cz < 9000 km s" 1 71 ±3 ± 7 

11 FP clusters, 1000 < cz < 11, 000 km s" 1 82 ±6 ± 9 

SBF for 6 clusters, 3800 < cz < 5800 km s" 1 70 ±5 ± 6 

4 Type II SN, 1900 < cz < 14, 200 km s" 1 72 ±9 ± 7 
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TABLE III Mean absolute B, V, and I magnitudes of nine SNe la with known Cepheid distances, without and with corrections for decline 
rate and color 



SN 


Galaxy 


(m — 


M)° 




My 


Mr 
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Ami 5 


(B - V)° 


M c R orr 


M£ orr 


M corr 
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(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


1937C 


1C 4182 
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-19.56 (15) 


19.54 (17) 
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-19.37 (17) 




1960F 


NGC 4496A 31.03 
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0.06 


-19.67 (18) 
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NGC 5253 
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Straight mean 
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Weighted mean 








-19.56 (07) 


-19.53 


(06) 


-19.25 (09) 








-19.47 (07) 


-19.46 (06) 


-19.19 (09) 



FIG. 1 Hubble Space Telescope in orbit photographed by the crew of the Space Shuttle Columbia, 
just after being released during mission STS-109 to service and upgrade Hubble. March 9, 2002. 



Credit: NASA, http://spaceflight.nasa.gov/gallery/images/shuttle/sts-109/html/sl09e5700.html 



FIG. 2 Comet P/Shoemaker-Levy 9 (1993e) "String of Pearls," HST/WFPC2, January 1994. 
Credit: NASA, H. A. Weaver (JHU) and T. E. Smith (STScI). 
http: / /hubblesite.org/ newscenter / archive /1994/13/ 



FIG. 3 Series of images showing the impact of Comet P /Shoemaker-Levy 9 (1993e) fragment W 
on Jupiter, HST/WFPC2, July 22, 1994. Credit: NASA and the HST Comet Impact Team. 



FIG. 4 The 
and 



Zahnlel p996) 



ight c urve from the impact of fragment R at 2.3 jx, adapted from lGraham et all (J1995) 



FIG. 5 The viewing geometry of a typical SL-9 impact (see text; adapted from IZahnld . 11996) 



FIG. 6 Flat Projection of Comet P/Shoemaker-Levy 9 (1993e) Fragments D (left) and G Impacts 
on Jupiter, HST/WFPC2, July 18, 1994. Two rings of propagating waves can be seen. Credit: 
NASA and the HST Comet Impact Team, http://hubblesite.org/newscenter/archive/1994/36/ 



FIG. 7 Globular Cluster 47 Tucanae (NGC 104), HST/WFPC2, July 1999. Credit: NASA and 
R. Gilliland (STScI). pittp://hubblesite.org/newscenter/archive/2000/33/ 



FIG. 8 Light curve of the slight dimming of the star HD 209458 due to a planet passing directly 
in front of it, HST/STIS, April-May 2000. Credit: NASA, T. M. Brown, D. Charbonneau, R. L. 
Gilliland, R. W. Noyes, & A. Burrows. 
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FIG. 9 Montage of four images showing jets from young stars HST/WFPC2. 
Top Left: HH 47 Credit: NASA and J. Morse (University of Colorado). 

http: / /hubblesite.org/newscenter/archive/199572471 To P Right: DG Tau B Credit: NASA, 



C. Burrows (STScI), J. Krist (STScI), K. Stapelfeldt (JPL), and the WFPC2 Science 
Team. http://hubblesite.org/newscenter/archive/1999/05/ Bottom Left: HH 34 Credit: 
NASA, J. Hester (Arizona State University), and the WFPC2 Investigation Definition 
Team. http://hubblesite.org/newscenter/archive/1995/24/ Bottom Right: HH 1-2 Credit: 
NASA, J. Hester (Arizona State University), and the WFPC2 Investigation Definition Team, 
http: / /hubblesite.org/ newscenter / archive / 1995 /24/ 



FIG. 10 HH 30 disk/jet, HST/WFPC2. Credit: NASA, C. Burrows (STScI), J. Krist (STScI), 
K. Stapelfeldt (JPL), and the WFPC2 Science Team, 
http: / /hubblesite.org/ newscenter / archive / 1999 /05/ 



FIG. 11 A bright one-sided microjet emerging from a proplyd located below Orion's Bright Bar. 
Adapted from Bally et al. (2000). 



FIG. 12 Protoplanetary disks (Proplyds) in the Orion Nebula (M42), HST/WFPC2. Credit: 
NASA, C. R. O'Dell (Vanderbilt University), and M. McCaughrean (Max-Planck-Institute for 
Astronomy), http://hubblesite.org/newscenter/archive/1995/45/ 



FIG. 13 Evolutionary tracks of 1 M Q , 5 M , and 25 M stars in the Luminosity-Effective Tem- 
perature (Hertzspring-Russell) diagra m. Thick segme nts mark long, nuclear burning, evolutionary 



phases (from Iben 1985; adapted from 



Prialnik 



2000) 



FIG. 14 The Luminous Blue Variable star Eta Carinae, HST/WFPC2, 1994. Credit: NASA and 
J. Hester (Arizona State University), http://hubblesite.org/newscenter/archive/1994/09/ 

FIG. 15 Rings around Supernova 1987A, HST/WPFC2, February 1994. Credit: NASA and C. Bur- 
rows (ESA/STScI). |http://hubblesite.org/newscenter/archive/1994/22/| 
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FIG. 16 Planetary Nebula MyCnl8, The Hourglass Nebula, HST/WFPC2. Credit: NASA, R. Sa- 
hai, J. Trauger (JPL) and the WFPC2 Science Team, 
http: / /hubblesite.org/ newscenter / archive / 1996 /07/ 



FIG. 17 Nebula surrounding the symbiotic star system He2— 104, the Southern Crab Nebula, 
HST/WFPC2, May 1999. Credit: NASA and R. Corradi (Instituto de Astrofisica de Canarias, 
Tenerife, Spain), M. Livio (Space Telescope Science Institute), U. Munari (Osservatorio As- 
tronomico di Padova-Asiago, Italy), H. Schwarz (Nordic Optical Telescope, Canarias, Spain), 
http: / /hubblesite.org/ newscenter / archive / 1999 /32/ 



FIG. 18 A schematic distribution of the planetary 



energy density (cr) plane (see text; adapted from 



nebulae morphology in the rotation (f2-magnetic 



Garcia-Seeura et al 



1999). 



FIG. 19 An HST-WFPC2 image of the circumstellar ring around supernova 1987A, obtained 
in May 2002 with narrow band filters that include Ha 6563 A and [Nil] 6584 A emission line 
radiation. We can identify a number of bright spots that correspond to gaseous protuberances 
on the inner side of the ring that are currently being hit by the supernova ejecta, thus becoming 
hot and luminous. Because of the narrow width of the filters employed to gather this image, the 
ring and the spots appear bright since they are mostly radiating in emission lines. On the other 
hand, both the supernova itself (a wide patch at the center of the ellipse) and a star (projected by 
chance near spot #6) are not visible in this image because most of their energy is radiated in the 
form of a rather smooth continuum. Credits: Nino Panagia (ESA/STScI), on behalf of the SINS 
Collaboration. 

FIG. 20 Jet and disk of stars and gas in the active galaxy M87, HST/WFPC2, 1994. Credit: 
NASA and H. Ford ( JHU) . |http : / /hubblesite.org/newscenter /archive / 1994/23/] 

FIG. 21 Black hole mass versus bulge lumin osity (left), and the luminosity-weighted dispersion 



(right). Adopted from 



Gebhardt et al 



(2000). 
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FIG. 22 Schematic gamma-ray burst from internal shocks and afterglow from external shock, 
arising from a relativistic jet produced by the collapse of a massive star. Internal shocks produce 
7-ra ys and neu t rinos, external shocks produce 7-rays, X-rays, optical and radio emission. Adapted 



from 



Msszamdtea). 



FIG. 23 Hubble Deep Field North, HST/WFPC2, December 1995. Credit: NASA, R. Williams 
(STScI) and the HDF Team, http://hubblesite.org/newscenter/archive/1996/01/ 

FIG. 24 Hubble Deep Field South, HST/WFPC2, October 1998. Credit: NASA, R. Williams 
(STScI) and the HDF Team. |http://hubblesite.org/newscenter/archive/1998/41/| 

FIG. 25 Optical and near- infrared images of distant galaxies in the Hubble Deep Field North. For 
each galaxy, the left panel shows a composite of WFPC2 images, while the right panel shows a 
new-infrared view with NICMOS. Courtesy of Mark Dickinson. 

FIG. 26 The star formation rate density versus r e dshift derived from the UV luminosity density. 



Adapted from 



Ferguson. Dickinson, and Williams! (J2000). 



FIG. 27 The redshift evolution of the co-moving stellar mass density. The vertical extent of the 
boxes shows the range of systematic uncertainty. The bottom two solid lines (on the right-hand 
scale) show the result of integrating the star-formation-rate histories traced by the rest-frame UV 
light, with and without corrections for dust extinction. The top two solid line s and the dashed 



line show theoretical predictions from semi-analytica 



Somerville. Primak. and Fabei 



2001). Adapted from 



galaxy evolution mod els ([Cole et al. 



Dicki nson et al. 



(2003). 



2000; 



FIG. 28 Joint confidence intervals for (Qm, ^a) from SNe la. Regions representing specific cos- 
mological scenarios are illustrated. 

FIG. 29 Hubble diagram of SNe la minus an empty (i.e., Qm + ^A = 0) universe comp ared to var 



ious cosmo 



( 1998) and 



ogical and astrophysics ! models. The points are redshift-binned data from 



Perlmutter et al 



Riess et al. 



( 1999). The observations of SN 1997ff are inconsistent wit h monotonic 



evolut ionary or dust effects that could mimic an accelerating universe. Adapted from 
( 20011 ). 



Riess et al. 
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